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Sound  Scattering  by  Cylinders 
of  Noncircular  Cross  Section 

Daniel  DiPerna 
Abstract 


This  thesis  investigates  a  new  method  for  predicting  the  farfield  scattered  pressure  of  a 
plane  wave  due  to  an  infinite  cylinder  of  noncircular  cross  section.  Both  impenetrable  and 
penetrable  fluid  boundary  conditions  will  be  treated  for  several  types  of  cross  sections 
and  for  a  large  frequency  range.  This  new  method  requires  the  conformal  mapping  of 
both  the  exterior  and  interior  of  a  closed  surface  to  a  semi-infinite  strip.  Numerically 
efficient  algorithms  will  be  presented  for  both  of  these  cases. 

A  new  method  for  satisfying  the  boundary  conditions  will  be  developed,  as  well 
as  an  efficient  algorithm  for  generating  the  required  modal  functions  on  the  boundary. 
Numerical  results  are  presented  for  cross  sections  in  the  shape  of  an  ellipse,  square,  and 
three  leaf  clover.  In  all  cases,  the  results  compare  extremely  well  with  exact  or  high 
frequency  asymptotic  results. 
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Chapter  1 


Introduction 


An  important  part  of  underwater  acoustics  research  deals  with  how  sound  interacts  with 
different  boundaries.  There  are  a  wide  variety  of  boundaries  that  exist  in  the  ocean 
environment.  There  are  two  dimensional  surfaces  such  as  the  ocean  surface  and  ocean 
bottom,  as  well  as  three  dimensional  boundaries  such  as  fish  or  submarines.  Another 
classification  of  the  boundary  is  the  material  composition  of  the  surface.  For  example, 
a  water  sediment  bottom  scatters  sound  differently  than  the  water  air  interface.  How 
sound  interacts  with  these  different  types  of  boundaries  is  important  in  the  design  of 
any  sonar  system.  It  is  important  to  quantify  how  these  different  boundaries  or  objects 
scatter  sound,  whether  it  be  for  the  detection  of  an  object  such  as  a  submarine,  or  for 
an  inverse  problem,  such  as  estimating  biomass  of  marine  organisms. 

The  manner  and  degree  to  which  an  object  scatters  sound  can  depend  strongly  on 
the  shape  of  the  scatterer,  as  well  as  on  the  relative  size  of  the  wavelength  of  the  sound 
wave  to  the  dimension  of  the  scatterer.  Exact  analytical  solutions  can  be  found  for  only  a 
limited  class  of  surfaces  where  separation  of  variables  is  applicable.  This  requires  that  the 
problem  be  formulated  in  a  coordinate  system  in  which  the  locus  of  points  corresponding 
to  the  radial  coordinate  equal  to  a  constant  coincides  exactly  with  the  scatterer  surface, 
and  that  the  Helmholtz  equation  in  that  coordinate  system  is  separable.  For  the  problem 


10 


treated  in  this  thesis-that  of  an  infinitely  long  cylinder  of  arbitrary  cross  section,  these 
requirements  are  met  only  for  three  cross  sections,  circular,  elliptic,  and  parabolic. 

If  separation  of  variables  is  not  applicable  for  a  particular  shape,  there  are  a  wide 
variety  of  methods  which  may  be  used.  All  of  these  methods  suffer  from  limitations  on 
either  the  frequency  range  or  the  class  of  surfaces  to  which  they  are  applicable.  For 
example,  one  analytic  approach  that  is  formally  exact  is  the  perturbation  method.  This 
may  be  used  for  penetrable  fluid  or  impenetrable  boundary  conditions,  but  is  only  valid 
if  the  shape  is  close  to  one  of  the  separable  geometries.  One  approximate  solution  is 
physical  optics.  However,  it  is  valid  only  at  high  frequencies  (i.e.  wavelength  much 
smaller  than  any  dimension  of  the  scatterer)  and  cannot  acount  for  the  penetrability  of 
the  scatterer. 

In  order  to  solve  a  problem  over  a  wide  range  of  frequencies  and  complex  shapes, 
one  must  resort  to  numerical  techniques.  Although  numerical  methods  which  are  for¬ 
mally  exact  have  been  developed,  they  can  suffer  from  computational  difficulties  when 
implemented.  For  example,  the  Transition  or  T-Matrix  method  is  a  numerical  technique 
which  is  used  extensively.  In  practice  when  the  surface  deviates  strongly  from  one  of 
the  separable  geometries,  the  resulting  matrix  that  needs  to  be  inverted  becomes  ill- 
conditioned.  Other  numerical  methods,  such  as  the  mode  matching  method  (MMM) 
may  be  numerically  intensive  for  even  intermediate  frequencies. 

As  can  be  seen  from  the  above,  there  are  many  different  approaches  applied  to 
the  problem  of  sound  scattering  by  cylinders  of  noncircular  cross  section.  In  general, 
problems  arise  when  the  scatterer  surface  diff''  significantly  from  a  geometry  in  which 
the  Helmholtz  equation  is  separable.  What  is  missing  is  a  numerically  efficient  method 
that  works  over  a  wide  frequency  range  for  a  large  class  of  surfaces  and  for  penetrable 
as  well  as  impenetrable  boundary  conditions.  The  method  presented  in  this  thesis,  the 
Fourier  Matching  Method  (FMM)  fills  this  void.  .An  overview  of  the  existing  methods 
and  the  FMM  is  given  in  the  following  table. 
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Matrix 

MMM 
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Near 
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All 
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All 

All 

Med. 

Med. 

Med. 

Restric- 

Near 

two- 
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None 

kS  «  1 

Separable 

None 

dimensional 

The  FMM  consists  of  the  following  procedure: 

1 .  Formulate  the  problem  in  a  new  coordinate  system  in  which  the  radial  coordinate 
u  =  constant  is  the  scatterer  surface, 

2.  express  the  solution  for  the  scattered  field  in  terms  of  cylindrical  eigenfunctions, 
and 

3.  solve  for  the  unknown  scattered  field  coefficients  by  expanding  the  boundary  con¬ 
ditions  in  a  Fourier  series,  and  then  forcing  the  first  M  terms  to  be  zero. 

Step  1  is  done  because  it  has  several  desirable  consequences.  By  using  a  special 
type  of  change  of  variables,  a  conformal  mapping,  calculation  of  the  normal  derivative 
of  the  pressure  (which  is  required  for  calculating  the  normal  particle  velocity)  becomes 
proportional  to  the  derivative  with  respect  to  u  of  the  pressure.  Another  benefit  is 
that  equispacing  sampling  points  in  the  new  angular  variable  v  samples  points  on  the 
surface  which  change  rapidly  at  a  higher  spatial  rate.  This  latter  benefit  has  the  effect  of 
requiring  fewer  sampling  points,  therefore  improving  the  numerical  efficiency.  The  only 
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drawback  to  the  conformal  mapping  is  that  since  conformal  mapping  is  a  two-dimensional 
procedure,  the  class  of  surfaces  that  the  FMM  can  treat  is  limited  to  two-dimensional 
scatterers. 

Step  2  assures  that  the  Helmholtz  equation  is  satisfied.  This  type  of  expansion  may 
be  used  independently  of  scatterer  shape,  as  opposed  to  the  T-Matrix  where  the  expansion 
functions  depend  on  shape.  The  benefit  of  expressing  the  solution  in  this  manner  is  that 
the  cylindrical  eigenfunctions  (Hankel  functions  and  complex  exponentials)  are  readily 
available  in  standard  math  libraries,  as  opposed  to  Mathieu  functions  (eigenfunctions  of 
the  Helmholtz  equation  in  elliptic  coordinates),  which  are  not  readily  available. 

Step  3  results  in  a  set  of  linear  equations  that  are  solved  for  the  unknown  coeffi¬ 
cients.  By  satisfying  the  boundary  conditions  using  this  manner,  the  number  of  coeffi¬ 
cients  required  to  achieve  a  converged  farfield  solution  is  much  smaller  than  the  MMM. 


1.1  Review  of  Relevant  Literature 

There  has  been  extensive  research  done  on  this  classic  problem  (for  example  see 
[1,  2,  3]  ).  This  section  will  give  a  brief  overview  of  the  existing  methods  that  deal  with 
the  problem  stated  above,  along  with  their  regions  of  applicability  and  shortcomings. 
The  three  distinct  types  of  methods  are  exact  analytic  solutions,  approximate  solutions, 
and  numerical  solutions. 


1.1.1  Analytic  Solutions 

The  scattering  of  sound  by  an  infinitely  long  cylinder  can  only  be  solved  exactly 
when  the  Helmholtz  equation  appropriate  for  that  geometry  is  separable.  This  means 
that  by  assuming  that  the  solution  is  a  product  of  a  function  of  the  angular  variable 
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and  a  function  of  the  radial  variable,  the  Helmholtz  equation  for  that  specific  geometry 
may  be  turned  into  two  ordinary  differential  equations.  Each  of  these  ordinary  differ- 
entiad  equations  depends  on  the  same  separation  constant,  and  may  be  solved  only  if 
the  constant  is  a  characteristic  value  or  eigenvalue.  There  are  infinitely  many  of  these 
eigenvalues,  with  each  eigenvalue  leading  to  a  solution  of  the  Helmholtz  equation.  One 
very  desirable  property  is  that  the  angular  eigenfunctions  form  a  complete,  orthogonal 
set,  so  any  function  may  be  represented  as  a  weighted  sum  of  these  functions.  The  total 
solution  then  is  an  infinite  summation  of  these  solutions,  with  an  expansion  coefficient 
2issociated  with  each  solution. 

Separation  of  variables  is  a  very  powerful  solution  technique,  but  in  two  dimensions 
the  number  of  separable  geometries  is  small,  namely  circular,  elliptical,  and  parabolic(4]. 
For  a  circular  cylinder,  the  problem  may  be  solved  exactly  for  impenetrable,  penetrable 
fluid,  and  elastic  boundary  conditions  [5,  6,  7,  8].  The  solution  of  the  ordinary  differential 
equations  in  this  coordinate  system  is  straightforward.  The  eigenvalues  are  simply  inte¬ 
gers,  leading  to  the  set  of  complex  exponentials  for  the  angular  eigenfunctions  and  Bessel 
functions  for  the  radial  eigenfunctions.  An  important  point  here  is  that  the  eigenvalues 
are  independent  of  w'avenumber,  which  means  that  even  for  penetrable  fluid  boundary 
conditons,  each  scattered  field  coefficient  is  coupled  only  to  one  incident  field  coefficient. 

For  elliptical  coordinates  there  is  only  a  one-to-one  coupling  between  scattered  and 
incident  field  coefficients  for  impenetrable  boundary  conditions.  Because  the  angular 
eigenfunctions  depend  on  wavenumber,  a  closed  form  expression  may  be  obtained  only 
for  impenetrable  scatterers  [9,  10].  For  penetrable  scatterers,  an  expression  involving 
the  inverse  of  a  matrix  has  been  obtained  by  Yeh  [11]  by  expanding  the  interior  field 
in  terms  of  the  exterior  eigenfunctions.  The  result  of  this  is  that  each  scattered  field 
coefficient  is  coupled  to  all  of  the  incident  field  coefficients.  Porter  [12]  has  obtained  an 
approximate  solution  for  weak  scatterers  with  small  eccentricity.  An  exact  solution  for 
elastic  boundary  conditions  has  not  been  found,  although  numerical  results  have  been 
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presented  by  Pillai  tt  al  [13]  using  the  T-Matrix  method.  In  Appendix  B,  separation  of 
variables  is  discussed  in  detail  for  circular  and  elliptic  coordinate  systems,  because  this 
method  provides  a  valuable  check  of  the  results  generated  by  the  method  presented  in 
this  thesis. 

Another  type  of  analytic  solution  is  a  perturbation  approach  [14,  15].  When  the 
scatterer  surface  deviates  only  slightly  from  circular,  the  boundary  conditions  may  be 
expanded  in  terms  of  a  Taylor  series.  By  equating  orders  of  the  perturbation  parameter 
<5,  correction  terms  to  the  zeroth  order  solution  (the  circular  solution)  are  obtained  to 
any  order  of  S.  As  in  the  case  of  the  penetrable  ellipse,  each  scattered  field  coefficient  is 
coupled  to  all  of  the  incident  field  coefficients,  even  for  impenetrable  scatterers.  These 
correction  terms  are  known  analytically.  The  drawback  of  this  method  is  that  the  per¬ 
turbation  series  will  only  converge  if  fC  1,  where  k  is  the  wavenumber  of  the  incident 
wave. 


1.1.2  Approximate  Solutions 

One  class  of  approximate  solutions  is  only  valid  at  high  frequency  [16].  One  such 
approximation  is  physical  optics[17,  18,  19].  This  approach  starts  with  a  Helmholtz 
integral  formulation,  and  assumes  that  each  point  on  the  surface  scatters  as  though  it 
were  part  of  a  plane  tangent  to  that  point.  The  resulting  integrals  may  be  evaluated 
asymptotically  by  stationary  phase  or  the  saddle  point  method.  This  is  a  very  powerful 
method  because  it  yields  insight  into  the  physics  involved  in  the  scattering  process, 
i.e.  it  shows  which  portions  of  the  object  contribute  strongly  to  the  scattered  field. 
The  limitations  are  that  it  cannot  account  for  the  penetrability  of  the  scatterer  or  for 
diffraction  from  a  sharp  edge. 

Another  high  frequency  approximation  which  may  take  into  account  sharp  edges  is 
the  geometric  theory  of  diffraction  [20,  21, 22,  23].  The  total  scattered  field  is  constructed 
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by  summing  contributions  from  speculary  reflected  rays,  surface  waves,  and  diffracted 
rays.  These  contributions  are  obtained  by  analyzing  canonical  cases.  For  example,  an 
edge  diffraction  contribution  is  calculated  by  analyzing  the  canonical  case  of  diffraction 
from  a  wedge. 


1.1.3  Numerical  Solutions 

There  are  various  numerical  solutions  for  the  scattering  of  sound  by  a  noncircular 
cylinder.  One  is  a  numerical  solution  of  the  integral  equations  [24,  25,  26].  An  integral 
equation  is  written  for  the  scattered  field  on  the  surface.  This  equation  is  then  solved 
by  discretizing  the  scattered  field  to  yield  a  system  of  linear  equations  that  may  be 
solved  for  the  unknown  surface  field.  The  major  drawback  of  this  method  is  that  for 
high  frequencies,  the  number  of  discretizations  and  thus  unknowns  becomes  very  large, 
making  the  method  computationally  intensive.  There  also  may  be  spurious  resonances 
corresponding  to  the  interior  eigenvalue  problem. 

Another  numerical  technique  is  the  extended  boundary  condition  or  T-Matrix 
method  [27,  3,  28,  13,  29].  This  method  assumes  an  expansion  of  the  surface  fields  in 
terms  of  a  complete  set  of  orthogonal  functions,  and  then  uses  the  Helmholtz  integral  for¬ 
mula  to  assure  that  the  Helmholtz  equation  is  satisfied.  The  result  is  that  the  boundary 
conditions  will  be  met  exactly,  but  in  truncating  a  matrix,  the  Helmholtz  equation  is  not 
satisfied  perfectly.  This  can  lead  to  problems  for  scatterers  whose  shapes  deviate  strongly 
from  circular,  as  the  matrices  that  need  to  be  inverted  tend  to  become  ill-conditioned. 
There  are  various  ways  to  deal  with  these  shortcomings.  For  example  Bates  and  Wall  [28] 
suggested  expanding  the  surface  fields  in  terms  of  a  more  appropriate  set  of  functions. 
For  an  elongated  scatterer,  it  is  desirable  to  expand  the  surface  fields  in  terms  of  elliptic 
cylinder  functions  (Mathieu  functions.)  Varadan  and  Varadan  [29]  have  also  suggested 
an  iterative  method  to  alleviate  the  ill-conditioning  for  elongated  scatterers. 


16 


Another  numerical  method  is  the  mode  matching  method  (MMM)  [30,  31].  The 
surface  fields  are  expanded  in  terms  of  a  series  of  solutions  (modal  functions)  of  the 
Helmholtz  equation,  with  an  unknown  coefficient  associated  with  each  modal  function. 
In  contrast  to  the  separation  of  variables  technique,  these  individual  solutions  are  not 
orthogonal.  This  makes  the  choice  of  the  unknown  scattered  field  coefficients  far  more 
difficult.  The  procedure  used  in  the  MMM  is  to  choose  the  unknown  coefficients  to  min¬ 
imize  the  mean  squared  error  in  the  boundary  condition.  This  minimization  is  achieved 
by  differentiating  the  mean  squared  error  on  the  surface  with  respect  to  each  scattered 
field  coefficient,  leaving  a  linear  system  of  equations  to  be  solved.  This  method  can  be 
proven  to  converge  to  the  correct  solution.  However,  a  large  number  of  terms  are  needed 
in  general  to  match  the  boundary  condition,  when  only  a  few  terms  are  necesary  to  cal¬ 
culate  the  farfield  scattered  pressure.  In  general,  most  of  the  computation  time  is  spent 
in  calculating  the  modal  functions,  so  there  is  an  inefficiency  in  carrying  more  terms  than 
are  necesary  to  calculate  the  farfield  scattered  pressure. 


1.2  Statement  of  the  Problem 

As  discussed  previously,  much  research  has  been  conducted  on  the  problem  of  sound 
scattering  by  an  infinitely  long  cylinder  of  noncircular  cross  section.  Each  of  these  is  valid 
for  a  very  specific  frequency  range,  boundary  condition,  or  class  of  surfaces.  In  this  thesis 
a  new  formulation  will  be  presented  that  will  efficient  ly  generate  farfield  scattered  pressure 
results  for  penetrable  fluid  or  impenetrable  infinitely  long  cylinders.  The  cross  sections 
of  the  cylinders  will  involve  ellipses,  a  square,  and  a  three-leaf  clover.  These  shapes  fall 
into  three  major  classes:  elongated  scatterers,  scatterers  with  sharp  edges,  and  scatterers 
with  concavities,  respectively.  The  formulation  is  the  same  regardless  of  surface  type, 
for  all  of  the  surfaces  mentioned  above,  with  the  surface  fields  being  expanded  in  terms 
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of  circular  eigenfunctions,  as  opposed  to  the  T- Matrix  method  where  the  choice  of  basis 
functions  depends  on  the  shape. 

The  fact  that  the  choice  of  basis  function  affects  the  numerical  properties  of  the 
T-Matrix  method  is  what  motivated  the  research  described  herein.  Bates  and  Wall 
suggested  that  the  ill- conditioning  problems  for  the  T-Matrix  could  be  alleviated  by 
choosing  a  coordinate  system  in  which  the  radial  coordinate  =  constant  is  as  close  to  the 
scatterer  as  possible,  and  expanding  the  surface  fields  and  Green’s  function  in  terms  of 
the  eigenfunctions  in  that  new  coordinate  system.  However,  this  still  assumes  that  the 
scatterer  surface  will  be  close  to  a  coordinate  system  that  is  separable,  and  requires  the 
calculation  of  special  functions  such  as  Mathieu  functions.  One  result  of  this  research 
is  that  by  expanding  the  surface  fields  in  a  Fourier  series  in  the  new  angular  variable  v 
and  the  Green’s  function  in  circular  eigenfunctions,  the  T-Matrix  will  be  able  to  handle 
even  elongated  scatterers  without  the  ill-conditioning  problems  usually  associated  with 
this  method. 

In  addition  to  using  the  circular  eigenfunctions  instead  of  more  complicated  special 
functions,  the  FMM  has  several  desirable  numerical  properties.  It  facilitates  and  makes 
more  efficient  the  calculation  of  the  circular  eigenfunctions  and  their  derivatives  on  the 
scatterer  surface.  It  also  makes  any  numerical  integration  over  the  scatterer  surface  more 
efficient  by  sampling  portions  of  the  surface  that  are  changing  rapidly  at  a  higher  spatial 
rate  than  portions  of  the  surface  that  are  smooth.  These  issues  are  discussed  in  Chapter 
2. 


The  major  contributions  of  this  thesis  are; 

1 .  An  efficient  algorithm  for  conformally  mapping  the  exterior  and  interior  of  an  ar¬ 
bitrary  shape  to  a  strip  is  developed  and 
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2.  a  new  method  for  satisfying  the  boundary  conditions  is  developed  which  improves 
the  numerical  efficiency  of  calculating  the  backscatter  from  noncircular  cylinders 
and  is  valid  for  a  large  class  of  cross  sections  and  frequencies. 


1.3  Outline  of  Thesis 

In  Chapter  2  the  class  of  surfaces  that  may  be  treated  by  the  method  presented  in 
this  thesis  is  discussed.  Algorithms  are  then  developed  to  conformally  map  the  exterior 
and  interior  of  a  closed  surface  to  a  semi-infinite  strip.  It  will  be  shown  how  the  Helmholtz 
equation  changes,  and  how  to  solve  this  new  Helmholtz  equation. 

In  Chapter  3,  various  numerical  issues  that  arise  in  the  solution  of  the  scattering 
problem  are  discussed.  These  issues  include  how  to  satisfy  the  boundary  condtions,  how 
to  truncate  the  resulting  infinite  series,  and  a  numerically  efficient  way  of  calculating  the 
modal  functions  and  their  normal  derivatives  on  the  surface. 

In  Chapter  4  the  scattering  results  generated  for  a  wide  variety  of  surfaces  and 
boundary  conditions  are  presented.  Included  are  results  for  high-aspect-ratio  elliptic 
cylinders,  a  square  cylinder,  and  a  cylinder  with  concavities.  Results  are  compared  to 
exact  solutions  where  applicable,  or  to  approximate  solutions. 

In  Chapter  5  there  is  a  discussion  of  the  method,  and  future  research  that  could 
follow  this  thesis  is  addressed. 
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Chapter  2 


Analytic  Development 


The  problem  treated  in  this  thesis  is  that  of  the  scattering  of  a  plane  wave  of  sound  by  a 
penetrable  fluid  cylinder  of  non-circular  cross  section.  The  acoustic  parameters  (density 
and  sound  speed)  of  the  obstacle  are  /Ji,Ci,  and  those  of  the  surrounding  fluid  are  p, c. 
The  geometry  is  shown  in  Fig.(2.1). 

The  approach  used  here  is  to  take  the  problem  formulated  in  Cartesian  coordinates, 
and  change  variables  (conformally  map)  to  a  new  coordinate  system  in  which  the  radial 
coordinate  =  constant  is  the  scatterer  surface.  This  has  the  effect  of  making  the  boundary 
conditions,  continuity  of  pressure  and  normal  particle  velocity,  easier  to  satisfy,  but  it 
also  makes  the  governing  equation  more  complicated.  However,  due  to  the  special  form 
of  the  change  of  variables,  the  solution  of  the  new  equation  is  also  obtainable. 


Consider  the  wave  equation 


1  d'^P 

c2 


(2.1) 


where  P  is  pressure,  is  the  Laplacian  operator,  c  is  the  speed  of  sound,  and  t  is 
time.  Using  the  form  of  the  Laplacian  operator  in  Cartesian  coordinates  and  assuming  a 
harmonic  time  dependence  where  u;  is  the  frequency  (and  suppressing  it  throughout 
this  thesis),  the  wave  equation  becomes  the  Helmholtz  equation: 


Pxx  "f  Pyy  +  P  =  0 


(2.2) 
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Figure  2.1:  Scattering  Geometry 


21 


where  k  —  uj/c  is  the  wavenumber.  Plane  waves  in  the  form  «»(<')+ v*in(M))  solutions 
of  this  equation,  where  is  the  angle  of  propagation  of  the  plane  wave  with  respect  to  the 
x-axis.  Each  plane  wave  is  associated  with  an  an  amplitude  A{fi),B{fj,).  By  integrating 
over  the  appropriate  contours  in  /i,  expressions  for  the  incident  and  scattered  field  may 
be  found. 

Changing  to  a  new  coordinate  system  changes  the  Helmholtz  equation,  giving  the 
following: 

Puu  +  Pvv  +  v)P  =  0  (2.3) 

where  (u,i;)  are  the  new  coordinates,  and  J-{u,v)  is  a  function  which  depends  on  the 
shape  of  the  scatterer.  Note  that  this  has  the  same  form  as  the  Helmholtz  equation 
in  Cartesian  coordinates  except  that  now  the  wavenumber  is  a  function  of  position. 
Therefore,  an  analogous  problem  is  wave  propagation  in  a  range  and  depth  dependent 
media.  This  is  a  very  difficult  problem  in  general  and  will  not  be  pursued  in  this  thesis. 
However,  by  solving  the  scattering  problem,  a  new  technique  for  solving  the  propagation 
problem  might  be  developed. 

Solutions  of  this  new  Helmholtz  equation  can  be  found  in  terms  of  plane  waves, 
but  now  the  wavenumber  is  a  function  of  position.  By  choosing  the  same  contours 
of  integration  as  in  the  circular  case,  solutions  may  be  found  in  terms  of  cylindrical 
eigenfunctions.  These  integral  representations,  coupled  with  the  fact  that  u  and  v  are 
orthogonal,  suggest  a  simple  and  efficient  w'ay  to  calculate  the  pressure  and  normal 
particle  velocity,  which  are  needed  to  satisfy  the  boundary  conditions. 

Changing  to  a  new  coordinate  system  also  generates  a  new  set  of  functions 
These  are  different  than  the  exponentials  that  would  ordinarily  be  generated  by  using 
cylindrical  coordinates  (e’”®).  The  new  angular  coordinate  v  is  a  more  natural  way  to 
describe  the  surface.  By  using  this  new  set  of  functions,  the  boundary  conditions  arc 
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satisfied  to  some  order  of  e*"",  as  well  as  satisfying  the  Helmholtz  equation  to  the  same 
order. 


The  remainder  of  this  chapter  is  divided  into  two  sections.  The  first  shows  how 
to  conformally  change  variables.  The  second  describes  the  resulting  Helmholtz  equation 
and  its  solutions. 


2.1  Conformal  Mapping 

In  this  section  it  is  shown  how  to  conformally  change  variables  to  a  new  coordinate 
system  in  which  the  radial  coordinate  =  constant  is  the  scatterer  surface.  This  has  several 
consequences.  One  is  that  the  conformal  nature  of  this  change  of  variables  ensures  that 
the  new  coordinates  are  orthogonal.  In  other  words,  by  changing  from  an  (x,  y)  coordinate 
system  to  a  (u,  v)  coordinate  system,  since  x  and  y  are  orthogonal,  u  and  v  will  also  be 
orthogonal.  This  is  a  desirable  property  because  it  facilitates  the  computation  of  the 
pressure  and  the  normal  particle  velocity  on  the  boundary,  which  are  needed  to  satisfy 
the  boundary  conditions.  Also,  by  expressing  the  field  quantities  in  terms  of  these  more 
natural  coordinates  instead  of  circular  cylindrical  coordinates,  points  along  the  surface 
which  change  rapidly  are  sampled  at  a  high  spatial  rate.  For  example,  for  a  square 
cylinder,  by  choosing  points  equally  spaced  in  angular  coordinate  v,  points  are  sampled 
more  rapidly  near  the  corners  where  the  field  is  expected  to  be  changing  more  rapidly. 

The  following  gives  a  procedure  for  generating  the  mapping  function  for  both  the 
exterior  and  the  interior  of  the  scatterer.  Once  these  functions  have  been  determined, 
the  uniqueness  of  the  map  is  checked  for  by  determining  where  the  Jacobian  of  the 
transformation  is  zero. 

Consider  the  class  of  surfaces  which  may  be  described  as  follows: 

^  =  Ps(<^)cos(</>) 
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y  =  Ps{<i>)  sin{4>) 


(2.4) 


where  p,,(f>  are  circular  cylindrical  coordinates.  The  only  constraint  is  that  ps{4>)  must 
be  single  valued,  i.e.  for  each  <^,  there  can  be  only  one  value  of  p,  (see  Fig. (2. 2)  for  an 
example  of  a  geometry  that  cannot  be  handled  by  this  conformal  mapping  algorithm). 

Without  loss  of  generality,  p,{4>)  may  be  expanded  in  a  Fourier  series  in  as  follows; 

OO 

p,{<f>)  =  a  +  ^  K  cos(n^)  +  r*  sin(n<?i.)]  (2.5) 

n=l 

where  a  is  the  average  radius,  and  and  r*  represent  the  deviation  of  the  surface  from 
circular.  Collecting  terms  this  may  be  rewritten  as 

OO 

X  +  tyUface  =  (2.6) 

n=l 

where 

=  5  K  +  ”-;i  (2.7) 

For  a  conformal  map  the  following  requirement  must  be  met: 

X  +  iy  =  M{z)  2  =  u  +  iv  (2.8) 

where  M(:)  is  an  analytic  function,  u  is  the  radial  variable,  and  i>  is  the  angular  variable. 
The  problem  therefore  is  to  choose  an  appropriate  M{z).  A  different  M{z)  is  needed 
for  interior  and  exterior  mappings.  These  mappings  are  different  because  one  of  the 
constraints  on  the  mapping  function  is  that  the  coordinate  system  becomes  circular  a-s 
the  radial  coordinate  increeises,  and  that  the  circles  must  be  outside  the  surface.  For  the 
interior  problem,  a  constraint  is  that  as  the  radial  coordinate  decreases,  the  coordinate 
system  becomes  circular,  and  that  the  circles  are  inside  the  surface.  Because  of  this 
difference,  there  are  different  forms  for  the  mapping  functions  for  the  interior  and  exterior 
problems. 
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2.1.1  Exterior  mapping 


Since  a  different  mapping  is  needed  for  the  interior  and  exterior  problems,  let  the 
mapping  function  M{z)  be  denoted  by  G(^)  for  the  exterior  problem.  The  function  G{z) 
must  be  chosen  so  that  the  following  three  conditions  are  met: 

1.  As  u  — ♦  oo,  the  coordinate  system  is  circular  (this  is  necesary  so  that  faxfield  results 
are  directly  comparable  to  existing  results  without  having  to  perform  a  possibly 
very  difficult  inverse  mapping), 

2.  the  transformed  Helmholtz  equation  is  solvable,  and 

3.  ti  =  0  is  the  scatterer  surface. 

Assume  the  following  form  for  the  mapping  function: 

OO 

G{z)  =  s-.e-  +  Y,  9ne-"-.  (2.9) 

n=0 

This  form  was  chosen  so  that  the  mapping  will  automatically  satisfy  Conditions  1  and  2. 
Condition  1  is  satisfied  because  as  u  — ^  oo,  e"”*  — >  0,  leaving  G{z)  —*  g-\e^  +  go,  which 
is  a  circular  coordinate  system  with  center  at  a:  =  Real{gQ),y  =  Imag{go).  Condition 
2  is  closely  tied  to  the  fact  that  the  mapping  must  be  unique,  which  will  be  discussed 
subsequently.  Condition  3  may  be  satisfied  bj'  choosing  the  appropriately.  To  do 
this,  the  original  parameterization  in  cylindrical  coordinates  of  the  surface  is  forced  to 
be  equal  to  the  mapping  function  when  the  radial  coordinate  is  zero,  or: 

OO 

X  +  iyUuTjace  =  G’(z)|„=o  =  5-1 (2.10) 

n=0 

Comparing  this  with  Eq.(2.6),  it  is  seen  that  there  is  a  fundamental  difference.  In 
Eq.(2.10)  there  are  only  negative  frequency  components  e'*”"  (other  than  e'"),  whereas  in 
Eq.(2.6)  there  are  both  positive  and  negative  frequency  components,  i.e. 
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From  this  it  is  concluded  that  (f>  ^  v.  Now  assume  that 


<f>{v)  =  u  +  sin(/t;)  +  cos(/t;)] .  (2.11) 

/=! 

This  form  was  chosen  so  that  if  the  surface  were  circular,  the  values  of  6  would  all  be 
zero,  and  (f>  =  v.  Since  the  problem  is  periodic,  we  expect  that  the  deviation  of  <t>  from  t; 
should  also  be  periodic,  and  therefore  be  representable  by  a  Fourier  series,  which  is  what 
is  inside  of  the  summation. 

The  problem  is  now  to  choose  the  values  of  S  so  that  Eq.(2.6)  and  E)q.(2.10)  are 
the  same.  Inserting  Eq.(2.11)  into  Eq.(2.10)  gives 


oo  oo 

+  Y.  (2.12) 


n=:l 


n=0 


Multiplying  both  sides  by  and  integrating  over  the  range  of  t;  leaves 


^  +  f]  1  dn  =  <  .  (2.13) 

JJ  [^,  j<l 

This  comes  from  the  orthogonality  relationships  of  the  complex  exponential  functions, 
i.e.  l/27r  =  6m,n,  where  is  the  Kronecker  delta  function.  The  upper 

result  is  a  set  of  non-linear  constraints  that  may  be  solved  for  the  values  of  S.  In  practice 
this  is  done  using  an  extension  of  the  Newton-Raphson  method  (Appendix  A).  Once  the 
values  of  6  are  known,  <?i(y)  is  also  known,  and  the  lower  result  is  used  to  calculate  gj. 

Once  the  mapping  coefficients  are  known,  the  uniqueness  of  the  map  must  be 
checked  for  (i.e.  that  each  (it,u)  corresponds  to  only  one  (x,  y)).  This  is  done  by  making 
sure  that  the  Jacobian  of  the  transformation  is  non-zero  in  the  desired  region  (outside 
the  scatterer.)  This  condition  is 

\G'{z)f^0,  u>0.  (2.14) 
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2.1.2  Interior  Mapping 


For  the  interior  mapping,  let  M{z)  be  denoted  by  F{z)^  and  z  =  ^  it},  where  ( 
is  the  radial  coordinate,  and  ij  is  the  angular  coordinate.  For  the  interior  mapping,  F(z) 
must  be  chosen  to  satisfy  the  following  conditions: 

1.  As  (  —*  — oo,  the  coordinate  system  is  circular, 

2.  the  transformed  Helmholtz  equation  is  solvable,  and 

3.  C  =  0  is  the  scatterer  surface. 

Now  assume  a  form  for  F(z)  as  follows; 

OO 

F(^)  =  E  /»«”'  (2-15) 

n=0 

Again,  Conditions  1  and  2  are  automatically  satisfied  by  the  form  of  F(z),  and  the  /„ 
must  be  chosen  to  satisfy  the  Condition  3.  Following  the  same  procedure  as  for  the 
exterior  problem,  let 

OO 

<f>(v)  =  ^  +  51  (f/  sin(/7;)  +  ef  cos(lTf)] ,  (2.16) 

/=i 

where  the  values  of  e  must  be  determined.  To  do  this  let 

OO  OO 

n=l  n=0 

Multiplying  both  sides  by  and  integrating  over  the  range  of  rf  leaves 

1  .jir  (  OO  0  7  >  1 

—  /  e‘^  \  ae‘^)  +  V]  \  drf  =  {  ~  .  (2.18) 

Jo  {  ^  I  fj  3  <  1 

Again,  the  values  of  c  are  solved  for,  and  then  values  of  fj  are  calculated  . 

For  this  map  to  be  unique,  the  Jacobian  must  be  non-zero  inside  the  scatterer,  or 
|F'(r)p  /  0,  C  <  0.  (2.19) 
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2.1.3  Example 


As  an  example,  let  p(^)  =  1  —  .2cos(2(^)  +  .1  sin(3(^)  4-  .1  cos(7^).  6j,  and  Cj  were 
calculated  for  values  of  j  =  40.  Figure(2.3a)  shows  lines  of  u  =  constant,  with  the 
innermost  curve  being  u  =  0,  along  with  the  parameterization  in  cylindrical  coordinates. 
As  u  increases,  the  coordinates  are  becoming  circular.  Figure(2.3b)  shows  the  surfaces 
C  =  constant,  with  the  outermost  curve  being  C  =  0,  again  with  the  originaJ  curve  in 
cylindrical  coordinates.  As  ^  decreases,  the  coordinates  are  becoming  circular.  Because 
the  coordinates  are  becoming  circular  far  away  from  the  boundary,  there  will  be  no 
need  to  perform  an  inverse  mapping.  As  can  be  seen  from  this  example,  there  is  excellent 
agreement  between  a  curve  given  in  cylindrical  coordinates  and  one  on  which  u  =  constant 
using  the  algorithm  given  above. 


2.2  Solutions  of  the  Helmholtz  Equation 

As  mentioned  previously,  changing  variables  conformally  changes  the  Helmholtz 
equation.  The  method  by  which  the  new  Helmholtz  equation  is  solved  is  described  in 
this  section.  The  method  is  broadly  similar  for  the  fields  outside  and  inside  the  scatterer. 
The  details  of  both  derivations  are  developed  in  the  following  subsections. 

2.2.1  Exterior 

After  conformally  mapping  the  exterior  of  the  scatterer  to  an  upper  halfspace,  the 
Helmholtz  equation  is  transformed  to  the  following; 

Puu  +  Pvv  +  k^Q{v.v)P  =  0  (2.20) 
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Figure  2.3;  Example  of  interior  and  exterior  mapping;  (a)  curves  u=  constant  for  several 
values  of  u,  with  the  inner  curve  showing  «  =  0  along  with  the  original  parameterization 
in  circular  coordinates,  (b)  curves  constant  for  several  values  of  u,  with  the  outer 
curve  showing  C  =  0  along  with  the  original  parameterization  in  circular  coordinates. 


where  G{u,v)  is  a  function  which  depends  on  the  specific  transformation.  More  specifi¬ 
cally,  it  is  the  Jacobian  of  the  transformation,  or 

Q{u,v)  =  xl  +  xl  (2.21) 

where  Xu  and  are  the  partial  derivatives  of  z(u,  v)  with  respect  to  u  and  t;.  Note  that 
the  uniqueness  condition  for  the  map  was  that  the  Jacobian  be  non-zero.  If  the  Jacobian 
were  zero,  this  would  correspond  to  a  zero  wavenumber,  or  an  infinite  soundspeed,  which 
is  non-physical. 


The  solutions  of  this  Helmholtz  equation  are  in  the  form 

gtA:[r  cos(/i)+v  sin(M)) 


(2.22) 


where  ^  is  a  parameter.  Physically,  this  solution  represents  a  wave  propagating  at  an 
angle  /x  with  respect  to  the  x  axis.  By  integrating  over  all  possible  values  of  //,  the 
solution  may  be  written  as 

P  =  f  (2.23) 

JWo  JW\ 

where  Ho  and  H^i  are  contours  chosen  such  that  the  integrals  will  represent  the  incident 
and  scattered  waves,  respectively,  and  A(fi)  and  B{n)  will  be  the  amplitudes  of  the 
incident  and  scattered  waves,respectively,  propagating  at  angle  fx.  In  general  A{fi)  is 
known,  and  B{fi)  is  a  function  that  must  be  solved  for.  Since  A{fi)  and  B{fi)  must  be 
periodic  functions,  they  may  be  expanded  as  Fourier  series: 

1  ^  ]  CJO 

n=— 00  n=:— 00 

Therefore,  the  values  of  a„  are  known,  and  the  values  of  are  unknown  constants 
describing  the  scattered  field. 


There  are  various  ways  to  evaluate  the  above  integrals,  each  with  different  numerical 
properties.  The  easiest  is  to  collect  terms  in  the  exponential  as  follows: 


r(u,i;)  =  y/x'^  +  y^ 

0(u,  u)  arctan  ^  .  (2.25) 
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Inserting  Eqs.(2.24  and  2.25)into  Eq.(2.23),  yields 


r  1  °° 

P  =  I  —  E 
Jwo  27r 


ginHgtfcr(u,i;)coa(6(u,w)— 


/■  1  °° 


(2.26) 

It  is  here  where  the  contours  necesary  for  Wq  and  Vl^i  become  apparent.  This  is  the 
same  integral  that  arises  in  the  solution  of  the  wave  equation  in  circular  coordinates. 
According  to  Sommerfeld  [32],  by  choosing  these  contours,  the  above  integral  may  be 
evaluated  analytically,  leading  to 


Where  and  are  Bessel  and  Hankel  functions  of  the  first  kind,  respectively,  of 
order  n.  Several  things  may  be  seen  immediately  from  this  expression.  If  the  scatterer 
is  circular,  0{u,v)  =  u,  and  r{u,v)  —  5-ie“,  and  the  solution  is  precisely  the  solution 
obtained  from  separation  of  variables  in  circular  coordinates.  For  the  non  circular  case, 
as  n  — >  oo,  d{u^v)  — >  u,  and  r(u,v)  —*  so  the  solution  will  approach  the  circular 

solution  (i.e.  a  summation  of  Bessel  and  Hankel  functions,  although  the  constants  b„ 
will  be  different)  far  from  the  scatterer.  The  first  term  in  this  expression  is  the  incident 
field,  and  the  second  is  the  scattered  field.  By  assuming  a  unit  amplitude  plane  wave 
excitation,  with  the  plane  wave  coming  from  an  angle  v,  with  respect  to  the  x-axis,  the 
total  field  \s  P  —  P"^^  +  , where  is  the  incident  field,  and  is  the  scattered 

field. 

OO 

=  J3  a„i"e‘”®<”’''’j„(itr(u,TO)  = 

n=  — OO 
OO 

pscat  ^  6„i"e’"^<’'’"’‘i/<*>(jfcr(n,u))  (2.28) 

n=— OO 


As  stated  above,  far  away  from  the  scatterer  the  coordinates  become  circular,  so 
using  this  fact  and  the  asymptotic  form  of  the  Hankel  functions,  the  farfield  scattered 


pressure  is  given  by 


pscat 
■*  u— *cx> 


(2.29) 


/“(v) 

where  f°°{v)  is  the  form  function[33].  This  is  the  quantity  that  will  be  generated  in  the 
results  section. 


2.2.2  Interior 


The  problem  of  determining  the  field  inside  the  cylinder  proceeds  in  almost  exactly 
the  same  manner  as  for  the  exterior  field.  Since  there  is  a  different  mapping  for  the 
interior,  there  is  also  have  a  different  Helmholtz  equation: 

Pi<  +  P,r,  +  k^^nC,rj)P  =  0  (2.30) 

where  ki  is  the  wavenumber  inside  the  scatterer,  and  P'{C,v)  is  the  Jacobian  of  the  inner 
mapping.  Following  the  same  procedure  as  for  the  exterior  case,  the  interior  pressure  is 

P=  f  (2.31) 

J  Wo  J  Wi 

where  C(//)  and  D{ii)  are  unknown  fields.  Expanding  each  of  these  as  Fourier  series  and 
integrating  leads  to 

oo  oo 

P=  E  Y.  (2.32) 

n=— OO  n=— oo 

and 


^(CiV)  Cretan  f— 'j 


(2.33) 


with  .V  and  y  expressed  in  terms  of  and  t}.  From  the  form  of  the  interior  mapping,  as 
C  — ♦  — 3C.  r((^,ri)  — >  0.  However,  the  Hankel  functions  become  infinite  as  its  argument 


34 


approaches  zero.  Since  the  field  must  be  finite  everywhere, 
interior  pressure 

i>"“=  ^ 

n=— oo 


=  0.  This  leaves  for  the 

(2.34) 


The  integral  which  led  to  the  above  expression  for  the  interior  pressure  is: 

pint  _  j  Y  (2.35) 

JWo  ntfc- 


Since  all  that  is  required  is  this  function  on  the  boundary,  let  [x{Co,Ti),y{^OiV)]  — 
[i(uo,t;),y(uo4v)].  Therefore,  the  interior  pressure  may  be  rewritten  as  a  function  of 
the  exterior  variable  u,v: 


n=— oo 


(2.36) 


where  p,  =  r(uo,  u),^^  =  u).  To  summarize,  in  this  section  expressions  were  de¬ 

veloped  for  the  pressure  fields  inside  and  outside  the  scatterer.  These  are  expressed  in 
terms  of  circular  eigenfunctions,  but  now  they  are  functions  of  the  new  angular  variable 
V.  In  Chapter  3  this  fact  will  be  exploited  in  developing  a  new  method  for  satisfying  the 
boundary  conditions. 


2.3  Boundary  conditions 

The  boundary  conditions  that  must  be  satisfied  on  the  scatterer  surface  are  con¬ 
tinuity  of  pressure  and  normal  particle  velocity.  In  the  remainder  of  this  section  all 
quantities  discussed  will  be  at  the  boundary.  Expressions  for  the  pressure  were  derived 
above.  The  normal  particle  velocity  V„  is  proportional  to  the  normal  derivative  of  the 
pressure  P„,  or 

K  =  ^  (2.37) 

iwp 

This  requires  the  calculation  of  the  normal  derivative  of  the  pressure.  Due  to  the  fact  that 
the  change  of  variables  was  conformal,  and  the  surface  is  u  =  Uo  ,  the  normal  derivative 
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is  easy  to  calculate.  The  normal  derivative  is  P„  =  n  •  VP,  where  n  =  is  a  unit  vector 
in  the  u  direction.  From  [34],  it  is  known  that 


du  hv  dv 


(2.38) 


where  d„  is  a  unit  vector  in  the  v  direction,  and  and  /i„  are  scale  factors  of  the 
transformation  and  are  given  by 


K  =  Jxl-k-  yl- 


(2.39) 


Therefore,  the  normal  particle  velocity  is  =  iwphu  •  Expressing  the  integrals  in  terms 
of  u  and  v  also  facilitates  the  computation  of  P„.  Consider  the  integral  representation 
of  the  typical  pressure  field  (this  is  for  the  scattered  field,  but  the  procedure  is  the  same 
for  the  incident  and  the  interior  field): 


1  .  oo 


Differentiating  with  respect  to  u  gives 


Rearranging  and  performing  the  integrations  yields 


where 


(2.40) 


=  ^  f  E  (x„  cos(/0  +  sin(,i))  (2.41 ) 

2n  Jw, 


OO 

Pu=  E  (2.42) 


^  =  2  f)  +  iyu(wo,  i^)) 

^  (.T„(«o.  v)  -  iy«(Mo,  t>))  • 


(2.43) 


So  now  expressions  have  been  derived  for  all  the  quantities  necesary  to  satisfy  the 
boundary  conditions,  namely  continuity  of  pressure  and  normal  particle  velocity.  The 
next  chapter  deals  with  the  numerical  issues  that  arise  when  actually  trying  to  satisfy 
the  boundai’}'  conditions. 
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Chapter  3 


Numerical  Considerations 


Some  of  the  numerical  issues  that  arise  in  the  solution  of  this  problem  are  treated  in 
this  chapter.  Among  these  are  how  the  boundary  conditions  are  satisfied,  a  convergence 
criterion,  and  how  to  calculate  the  basis  functions  efficiently. 


3.1  Boundary  Conditions 

The  most  difficult  part  of  the  scattering  problem  under  consideration  is  how  to 
satisf}'  the  boundary  conditions.  Although  the  expressions  that  were  derived  for  the 
scattered  field  on  the  surface  appear  to  be  similar  to  that  of  the  circular  cylinder,  there 
is  one  major  difference:  the  modal  functions  are  not  orthogonal.  This  makes  it  much 
more  difficult  to  choose  the  scatered  field  coefficients  properly.  In  the  following  section 
the  problem  will  be  explained  better,  and  it  will  be  shown  how  to  overcome  it. 

For  simplicity,  the  impenetrable  case  will  be  discussed,  but  the  issue  is  the  same 
for  the  penetrable  case.  The  expressions  for  the  fields  derived  in  the  previous  chapter 
are  expanded  in  terms  of  basis  functions  which  are  solutions  of  the  Helmholtz  equation. 
For  example,  from  Eq.(2.28)  the  basis  functions  Pm{i')  for  the  scattered  pressure  on  the 


37 


surface  are 


(3.1) 


When  treating  rigid  boundary  conditions,  the  following  b<isis  functions  from  Eq.(2.37 
and  2.42)  will  be  used: 


1  dPm 
iujphu  du 


(3.2) 


An  infinite  summation  of  these  basis  functions  should  be  able  to  match  the  bound¬ 
ary  condition.  In  practice  however,  infinite  summations  cannot  be  performed,  and  there¬ 
fore  the  summations  must  be  truncated  at  some  point.  This  is  a  major  obstacle,  because 
it  is  far  from  obvious  where  to  truncate  this  summation.  For  the  separable  geometries 
the  modal  functions  are  eigenfunctions  of  a  Sturm- Liouville  problem,  and  therefore  are 
orthogonal.  Multiplying  the  boundary  condition  equation  by  one  of  the  eigenfunctions 
and  integrating  over  the  angular  variable  will  thus  leave  a  closed  form  expression  for  the 
unknown  constant.  Therefore,  truncating  the  summations  for  a  separable  geometry  will 
not  affect  the  accuracy  of  the  constants  that  are  solved  for. 


For  a  nonseparable  geometry,  the  modal  functions  are  no  longer  orthogonal.  Be¬ 
cause  of  this  fact,  questions  arise  on  which  functions  are  appropriate  by  which  to  multiply 
the  boundary  condition  (since  there  are  no  eigenfunctions).  No  matter  what  set  of  func¬ 
tions  are  chosen,  the  result  of  multiplying  by  the  weighting  functions  and  integrating 
will  be  a  system  of  equations  to  be  solved  for  the  unknown  coefficients.  It  is  here  where 
the  truncation  problem  arises.  Since  all  of  the  coefficients  are  coupled  to  each  other, 
truncation  of  the  series  will  affect  the  accuracy  of  the  constants  that  are  solved  for. 

The  following  example  will  demonstrate  how  this  non-orthogonality  of  modal  func¬ 
tions  can  make  the  boundary  conditions  far  more  difficult  to  satisfy. 
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3.1.1  Example:  Circular  cylinder 


Before  treating  the  noncircular  cylinder,  a  case  will  be  examined  for  which  an  exact 
solution  exists,  the  circular  cylinder.  In  this  case,  separation  of  variables  shows  that  the 
basis  functions  on  the  surface  of  the  scatterer  are 

VCM  =  (3.3) 

tupn^ 

This  is  the  standard  expression  since  for  a  circular  coordinate  system  v  =  6  and  = 

a,  where  a  is  the  radius  of  the  scatterer.  Here  the  basis  functions  are  proportional  to 
which  constitute  an  orthogonal  set,  and  therefore  any  incident  field  may  be  represented 
in  terms  of  this  set  of  basis  functions.  The  truncation  properties  of  this  problem  are  well 
known.  The  series  may  be  truncated  so  that  (m(  <  A/,  where  M  is  somewhat  larger 
than  ke^°  .  This  truncation  point  is  arrived  at  because  the  incident  field  is  expanded  in 
terms  of  J^{he'^).  Once  M  >  A:e“°,  the  value  of  the  bessel  function  becomes  small,  and 
therefore  the  amount  of  energy  contained  in  this  mode  is  negligible.  Therefore,  in  terms 
of  a  finite  summation,  the  boundary  condition  becomes 

M 

L  +  =  (3.4) 

m=— A/ 


where  is  the  known  incident  particle  velocity  on  the  boundary  and  the  6^  are  the 

scattered  field  coefficients  that  must  be  chosen.  This  would  be  solved  by  multiplying  by 
g-.mv  ^  integrating  over  the  range  of  c  to  get 


hm  —  Q 


where  =  e  and  Oq  is  the  angle  of  incidence  of  the  incident  plane  wave. 


Now  suppose  that  the  coefficients  are  a  linear  combination  of  some  new  set  of 
coefficients 

b,r,  =  f;  (3.()) 
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where  Q  is  some  known  matrix.  Inserting  this  into  Eq.(3.4)  and  rearranging  the  summa¬ 
tions,  leaves 

1"'“  =  -  Y.  ‘S.*’  E  •  (3.7) 

n=—N  m=—M 

' - - - ' 

This  has  the  effect  of  expanding  the  solution  in  terms  of  a  new  set  of  basis  functions 
Note  that  these  functions  are  no  longer  orthogonal,  although  they  are  each  still 
a  solution  of  the  Helmholtz  equation,  because  each  is  a  linear  combination  of  the 

Further  assume  that  N  is  less  than  M.  Under  this  assumption,  trying  to  satisfy 
the  boundary  conditions  in  the  same  manner  as  before  results  in  an  overdetermined 
problem.  Multiplying  by  e"”"’'  for  |m|  <  M  and  integrating  gives  M  constraints,  amd 
only  N  unknowns.  This  may  be  solved  by  minimizing  the  mean  squared  error,  but 
nonetheless,  error  will  be  introduced.  This  method  will  be  referred  to  as  Method  1. 

Another  alternative  is  to  multiply  by  for  \m\  <  which  would  yield  an 

exactly  determined  problem,  but  now  there  is  no  constraint  on  the  c’""'  for  |A^|  <  m  <  |Af  | 
so  again  there  will  be  error.  This  will  be  referred  to  as  Method  2.  The  conclusion  therefore 
is  that  in  either  case  there  will  be  some  error  introduced,  although  the  error  is  not  the 
same  in  each  case.  In  Method  1  the  error  will  be  distributed  throughout  all  of  the 
angular  frequency  components  e'"*",  whereas  in  Method  2  all  the  error  will  be  at  the 
higher  frequencies. 

From  this  example  it  is  seen  that  by  choosing  the  incorrect  set  of  basis  functions  a 
problem  that  can  be  solved  exactly  can  be  turned  into  one  that  will  necesarily  contain 
some  error. 

It  should  not  be  a  surprise  that  a  problem  that  can  be  solved  exactly  can  be  turned 
into  a  more  difficut  one.  In  this  case,  choosing  a  specific  Q  {Qn,m  = 

physically  represents  moving  the  coordinate  origin  from  the  center  of  the  circle  to  Xc  — 
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Ti  cos(t;i ),  xfc  =  ri  sin(t;j).  The  result  of  this  is  a  radius  that  changes  as  a  function  of  angle, 
which  is  the  same  as  the  noncirculax  case.  Therefore,  it  is  expected  that  the  resulting 
equations  should  be  similar  to  that  of  the  noncircular  case. 


3.1.2  Noncircular  cylinder 

The  same  problem  of  non-orthogonality  of  modal  functions  arises  in  conjunction 
with  the  noncircular  cross  section  problem.  There  is  an  infinite  summation  of  non- 
orthogonal  basis  functions,  which  must  be  truncated  somewhere.  However,  once  the  series 
is  truncated  error  will  be  introduced  when  choosing  the  unknown  coefficients,  presumably 
because  an  inappropriate  set  of  basis  functions  was  used.  Ideally,  there  would  be  some 
way  to  construct  the  appropriate  functions  from  the  basis  functions  derived  previously  in 
conjunction  with  the  transformed  Helmholtz  equation,  but  this  is  a  very  difficult  problem. 
Since  it  is  not  known  a  priori  what  the  appropriate  basis  functions  are,  a  method  must 
be  developed  to  choose  the  unknown  coefficients  appropriately. 

Both  methods  1  and  2  described  above  are  in  fact  special  cases  of  the  same  method. 
Recall  the  boundary  condition  after  truncating  the  series  for  the  scattered  field  is 

M 

£  +  »''“(>')  =  0-  (3.8) 

m=-M 

Multiplying  by  a  set  of  functions  and  integrating  over  the  range  of  v  yields 

f2ir  /•2r 

Y.  Vr\v)Mv)dv+  V"»0„(T;)du  =  0  (3.9) 

m=-M 

for  each  n.  The  question  now  is  to  choose  an  appropriate  set  of  xj).  A  logical  choice  is  to 
try  to  minimize  the  mean  squared  error  fl  on  the  surface,  or 

II  ^  II 

n=  E  (3.10) 

|m=— A/  II 
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where 


M 

M 

2 

L 

dv 

m=—M 

ms— M 

E  +  =  /  E  +  dv  (3.11) 

||m=— Af  II  ®  m=— M 

In  order  to  minimize  this  expression,  differentiate  with  respect  to  6„  and  set  equal  to  zero 
to  obtain 

r  E  b^Vr\v)Vr'-iv)dv+  rV'-^{v)Vr^-iv)dv=  (3.12) 

for  n  <  M.  Comparing  Eq.(3.12)  to  Ek}.(3.9),  it  is  seen  that  this  is  the  same  as  choosing 
=  V^*“‘*(n).  This  choice  of  weighting  function  is  equivalent  to  the  mode  matching 
method  (MMM)[31].  There  is  a  proof  that  this  will  converge  to  the  true  solution  with 
increasing  M,  but  in  practice,  M  will  be  very  large. 


Another  choice  for  ^  is  to  let  =  e*”".  This  choice  of  weighting  function  is  referred 
to  as  the  Fourier  matching  method  (FMM).  It  is  here  where  the  conformal  mapping 
becomes  important.  Recall  that  this  v  is  specific  to  the  geometry.  It  w'as  generated 
through  the  conformal  mapping  algorithm.  Without  changing  variables  to  (u,  v),  this 
choice  of  xj)  would  not  be  possible.  The  author  cannot  prove  rigorously  that  this  will 
converge  to  the  true  solution.  However,  as  can  be  seen  in  Chapter  4,  in  all  of  the  ceises 
that  have  been  tried,  the  farfield  scattered  pressure  calculated  using  the  FMM  compares 
extremely  well  with  either  the  exact  solution  (for  the  elliptic  cylinder),  the  high  frequency 
limit  predicted  using  the  physical  optics  approximation  (for  the  impenetetrable  ellipse 
and  three-leaf  clover),  or  the  DWBA  (for  the  weakly  scattering  square).  The  FMM  is  in 
general  preferable  to  the  MMM  for  calculating  the  farfield  scattered  pressure  because  far 
fewer  terms  are  required  for  a  convergent  solution.  For  example,  for  a  2:1  aspect  ratio 
rigid  elliptic  cylinder  using  the  MMM,  M  =  50  is  required  for  a  convergent  solution  at 
ka  =  3.  For  the  same  problem  using  the  FMM,  there  is  a  converged  solution  for  M  =  12. 
However,  if  a  nearfield  solution  is  desired,  the  MMM  must  be  used,  since  the  error  in 
satisfying  the  boundary  condition  using  the  FMM  does  not  approach  zero  with  increasing 


In  summaxy,  the  equations  for  the  boundary  conditions  using  the  FMM  are  for  rigid 
boundary  conditions: 

r2ir  M  ,2)r 

E  ' 

m=—M 

for  soft  boundary  conditions, 


r  E  =  ^  (3.13) 

•'0  Jo 


6„P'“‘(t;)c-’‘”do  +  r  =  0,  (3.14) 

Jo 


/.2»  M 
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m=—M 

and  for  penetrable  fluid  boundary  conditions, 


f7ir  M  f2ir  .  .  flu  M 

/  E  /  P‘"'(t;)c-"Vt;  =  /  ^  Cr^P'^\v)e-'’"'dv 


2ir  M 

E 

m=—M 


fir  "U  r2ir  -  fir 

/  E  /  V‘’*»e— dt;  =  ^/  c„.V-‘(t;)e-‘-dt; 

*'0  m=-M  m=-A/ 

(3.15) 


for  n  <  M. 


3.2  Convergence 

Since  for  the  FMM,  the  error  on  the  boundary  does  not  necesarily  decrease  with 
increasing  M,  another  criterion  for  convergence  is  used.  Therefore,  a  criterion  using  the 
field  quantities  away  from  the  scatterer  surface  will  be  used.  The  following  is  a  measure 
of  how  well  the  optical  theorem  is  satisfied  [31],  and  is  used  in  the  MMM.  Since  the 
MMM  and  the  FMM  have  similar  forms,  it  may  be  used  directly  term  by  term  as  given 
in  [31]: 


e(M)  =  2 


(3.16) 


|E;if=-A/  K?  -  ReYJ:i=.M 

where  ({M)  is  the  error  as  a  function  of  the  number  of  unknown  coefficients.  The  nu¬ 
merator  is  a  measure  of  how  well  the  optical  theorem  is  satisfied.  The  first  term  in  the 
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numerator  is  simply  the  total  cross  section  a,  which  is  the  total  scattered  power  normal¬ 
ized  by  the  incident  power  density.  The  denominator  is  approximately  twice  the  total 
cross  section  a  and  therefore  normalizes  the  error  in  the  optical  theorem  by  the  totail 
scattered  power. 

In  practice,  M  will  be  increased  until  e(M)  is  smaller  than  some  given  tolerance. 
This  is  a  necesary  condition  that  the  field  must  satisfy  (i.e.  this  is  a  statement  of  con¬ 
servation  of  energy,  so  it  must  be  satisfied),  but  it  is  not  suflScient. 

The  insufficiency  of  this  criterion  occurs  because  there  are  many  choices  for  6^  that 
satisfy  the  above  condition.  If  the  6^  were  to  be  chosen  simply  to  satisfy  this  condition, 
the  choice  would  not  be  unique.  The  above  condition  is  totally  independent  of  the  type 
of  scatterer,  both  shape  and  whether  the  boundary  condition  is  soft,  rigid  or  penetrable. 

The  farfield  condition  becomes  useful  only  when  used  in  conjunction  with  the  con¬ 
ditions  imposed  by  the  FMM  on  the  scatterer  surface.  As  stated  previously,  the  FMM 
only  matches  the  boundary  condition  to  some  order  of  e*”*".  It  is  unclear  how  to  choose 
M.  The  above  criterion  gives  us  a  measure  of  how  large  M  should  be.  For  example,  the 
boundary  condition  may  be  met  in  the  FMM  sense  for  some  M.  If  that  M  is  not  large 
enough,  may  not  be  small,  and  therefore  that  choice  of  will  not  satisfy  a  necesary 
condition.  By  increasing  M  until  cm  is  less  than  some  tolerance.  Therefore,  out  of  all 
the  possible  choices  of  bm  that  satisfy  the  above  criterion,  we  choose  the  ones  that  also 
satisfy  the  boundary  conditions  in  the  FMM  sense. 

In  all  caises  that  have  been  tried,  when  c(A/)  decreases  to  within  some  suitably  small 
tolerance,  the  backscattered  pressure  is  no  longer  changing  significantly,  and  compares 
well  with  the  reference  solutions.  Specific  examples  are  given  in  Chapter  4. 
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3.3  Numerical  Computation  of  Modal  Functions 


In  Chapter  2,  expressions  were  derived  for  the  pressure  and  velocity  fields  in  terms 
of  series  of  Bessel  functions.  When  trying  to  satisfy  the  boundary  conditions,  2M  +  1 
Bessel  functions  must  be  calculated  for  N  integration  points.  This  can  be  numerically 
intensive.  The  following  gives  an  improved  method  for  calculating  the  mode  functions. 
Recalling  the  integral  representation  of  the  scattered  field  and  rearranging  somewhat 


leaves 


where 


pscat^  I  g 
- - 


A'(u,u,/x)  = 


(3.17) 


(3.18) 


Since  /\(u,t;,/i)  is  a  periodic  function  in  fi,  it  may  be  expanded  in  a  Fourier  series  as 


follows: 


1  °° 

K(u,v,n)  =  — 

H  — —  OO 


where  k„  is  a  function  of  u  and  v,  and  may  be  found  by  integrating  as  follows 


Jo 


(3.19) 


(3.20) 


The  Kn  are  known  analytically,  but  in  genera]  the  result  is  a  multiple  summation  of 
products  of  Bessel  functions,  and  the  truncation  of  these  series  is  far  from  obvious.  It 
is  also  far  more  efficient  to  perform  the  above  integration  numerically,  and  since  that 
integration  is  essentially  a  Fourier  transform,  a  Fast  Fourier  Transform  (FFT)  algorithm 
may  be  used  to  calculate  it.  For  simple  cases,  where  there  is  only  one  roughness  coefficient, 
it  may  be  easier  to  find  «„  analytically.  For  example,  for  an  ellipse,  gi  is  the  only  nonzero 


coefficient,  and 


«n  =  J„{kg,e-nre'^^. 


(3.21) 


In  general,  the  resulting  integral  becomes 


f  [  1  r  ^  1 

pscat  ^  I  Y  (3.22) 

Ln=— oo  171  =  — X 
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Multiplying  these  series  together  gives 


OO  OO  f 

Y.  Y  (3.23) 

n=-oo/=-oo 

This  may  be  evaluated  analytically,  and  after  rearranging  the  summations, 

E  E  « _ e").  (3.24) 

m=~oo  n=:~oo 

Although  this  appears  more  complicated  than  the  previous  expression  (Eq.(2.28)),  it  is 
more  efficient  to  evaluate  because  each  needs  to  be  evaluated  only  once  for  a  given 
frequency,  instead  of  N  times  for  the  previous  expression. 


A  similar  expression  may  be  derived  for  the  normal  particle  velocity.  Due  to  the 
conformal  nature  of  the  change  of  variables,  ^  Therefore  the  normal  particle 

velocity  is  Performing  this  differentiation  yields 

yacat  _ — f  ^  6„e‘"'‘e‘^®~*'““*^*'~''^A'(ii.v,//)(x„cos(//) +  i/„sin(/i)]d/i.  (3.25) 

txjph^  JWi 

By  letting 

A"(u,  v,p)  =  K (u,  V,  fi)  [x„  cos(/i)  +  Vu  sin(M)] ,  (3-26) 


K'  can  be  expanded  as  a  Fourier  series: 

CC 

K'(u,v,^)=  Z  (3-27) 

n=— OO 

Following  the  same  procedure  as  for  the  pressure  the  normal  particle  velocity  is  given  by: 

(3.28) 

The  normal  particle  velocity  of  the  incident  field  is  known  immediately  from  Eq.(2.37): 

(3.29) 


1  OO  OO 

E  6„  E 

m=—oo  n=— OO 


[x„  cos(n,)  +  yusin{vi)] . 

uphu 


The  integral  which  led  to  the  above  expression  for  the  interior  pressure  is 


pxnt 


(3.30) 
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Following  the  same  procedure  as  for  the  exterior  problem,  gives 


where 


n=-oo 


T(u,t;,p)  =  e**‘l*co.(M)+y.m(M)l 


(3.31) 


(3.32) 


Since  T{u,v,fi)  is  a  periodic  function  in  /it,  it  may  be  expanded  in  a  Fourier  series  as 


follows: 


r(x,  ^  E 

n=-“Oo 


where  is  a  function  of  u  and  w,  and  may  be  found  by  integrating  as  follows 


rn  =  /  e  ‘"'‘T(tt,t;,/t)d/t 
Jo 


(3.33) 


(3.34) 


Inserting  Eq.(3.34)  into  Eq.(3.31)  in  and  multiplying  the  series  yields 


oo  oo  . 

E  E  e‘"-d/x. 

n^r'oom^oo 


(3.35) 


This  integral  is  zero  unless  n  =  0,  so  letting  n  =  0  and  rearranging  the  sum,  the  internal 
pressure  becomes 

OO 

E  (3.36) 

m=— OO 

The  normal  particle  velocity  may  also  be  found  very  easily.  Taking  the  normal 
derivative  of  the  internal  pressure  with  respect  to  the  normal  gives: 


I  r  °° 

=  —  /  y 

^phu  Jwo 


c^e'''^T'{u,v,n)dft 


where 


T'{u,v,n)  =  T{u,v,fi)  [x„cos(/t)  +  y„sin(/t)] 


(3.37) 


(3.38) 


T'  can  be  e.xpanded  as  a  Fourier  series: 


T\xi,v,n)= 


(3.-39) 
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Inserting  Eq.(3.39)  into  Eq.(3.37)  and  integrating,  yields 

E  (3.40) 

oo 

A  numerically  efficient  way  to  generate  all  of  the  quantities  that  are  needed  to 
satisfy  penetrable  fluid  boundary  conditions  at  the  surface  of  the  scatterer,  the  interior 
and  exterior  pressure  fields,  and  the  interior  and  exterior  normal  particle  velocities  has 
been  developed.  This  representation  is  numerically  efficient  for  several  reasons.  One 
is  that  Hankel  functions  no  longer  have  to  be  calculated  for  every  integration  fKjint 
along  the  boundary.  A  second  is  that  by  equally  spacing  integration  points  in  the  new 
angular  variable  u,  points  along  the  surface  that  change  rapidly  are  sampled  at  a  higher 
spatial  rate,  so  fewer  integration  points  are  required  than  if  points  were  equally  spaced 
in  cylindrical  coordinates  . 

It  has  also  been  discussed  how  to  choose  the  scattered  field  coefficients,  and  a 
convergence  criteria  for  the  farfield  scattered  pressure. 
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Chapter  4 


Numerical  Results 


This  chapter  will  begin  by  comparing  the  numerical  efficiency  of  the  FMM  with  various 
existing  methods.  Numerical  results  generated  using  the  FMM  will  then  be  presented 
for  a  several  types  of  cross  sections  and  boundary  conditions  and  compared  with  exact 
solutions,  where  applicable,  or  approximate  solutions.  All  results  involve  plots  of  the 
backscattered  form  function  as  defined  in  Eq.  (2.29)  versus  a  dimensionless  term  k 
times  a  characteristic  dimension  of  the  cross  section.  The  equations  used  to  generate  all 
results  for  the  FMM  are  given  by  Eq.(3.13)  for  rigid  boundary  conditions,  Eq.(3.14)  for 
soft  boundary  conditions,  and  Eq.(3.15)  for  penetrable  fluid  boundary  conditions.  The 
modal  functions  for  the  scattered  and  interior  fields  are  those  given  in  Eqs.(3.2 '  3.28) 
and  Eqs.(3.36  and  3.40),  respectively. 


4.1  Numerical  Efficiency 

Two  numerical  methods  that  are  used  to  solve  the  problem  of  scattering  by  a 
noncircular  cylinder  are  the  MMM  and  the  T-Matrix  method.  One  similarity  between 
both  of  these  approaches  and  the  FMM  is  that  they  all  require  the  calculation  a  Hankel 
function  (for  each  unknown  coefficient)  at  many  points  along  the  boundary,  and  then 
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an  integration  over  the  boundary.  Therefore,  the  numerical  efficiency  of  each  of  these 
methods  depends  on  two  criteria: 

1 .  The  number  of  unknown  coefficients,  and 

2.  the  number  of  integration  points. 

For  smooth  scatterers  that  do  not  deviate  strongly  from  circular,  there  will  not  be 
a  significant  difference  between  the  numerical  efficiency  of  the  FMM  and  the  T- Matrix  or 
the  MMM.  However,  in  more  extreme  cases,  the  FMM  can  be  much  more  efficient  than 
the  other  approaches. 

For  example,  the  FMM  is  numerically  more  efficient  than  the  T- Matrix  for  elon¬ 
gated  cross  sections  because  fewer  terms  are  needed  for  the  numerical  integrations.  This 
arises  due  to  the  fact  that  the  conformal  mapping  used  in  the  FMM  samples  points  along 
the  scatterer  surface  that  change  rapidly  at  a  higher  spatial  rate.  The  following  table 
shows  the  number  of  modes  and  integration  points  required  for  the  FMM  and  the  T- 
Matrix  for  the  c<ise  of  a  rigid  10:1  aspect  ratio  elliptic  cylinder  for  kL  =  5,  where  L  is 
the  length  of  the  semimajor  axis.  Note  that  in  order  to  make  the  T-Matrix  stable,  the 
basis  functions  for  the  surface  field  were  expanded  in  terms  of  a  Fourier  series  in  u,  and 
the  Green’s  function  in  cylindrical  eigenfunctions.  The  MMM  required  more  terms  than 
the  computer  could  handle: 


Method 

FMM 

T-Matrix 

MMM 

Number  of  Modes 

18 

18 

»  18 

Number  of  Integration  Points 

120 

800 

800 

The  FMM  is  numerically  more  efficient  than  the  MMM  because  far  fewer  unknown 
coefficients  are  required  for  the  same  degree  of  accuracy.  The  MMM  requires  more  terms 
is  because  it  strictly  matches  the  boundary  condition,  whereas  the  FMM  may  only  be 
used  to  calculate  the  far  field  pressure.  For  example,  the  following  table  shows  the  number 
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of  modes  and  integration  points  required  for  the  FMM,  the  T-Matrix  and  the  MMM  for 
the  case  of  a  rigid  2:1  aspect  ratio  elliptic  cylinder  for  kL  =  3: 


Method 

FMM 

T-Matrix 

MMM 

Number  of  Modes 

12 

12 

52 

Number  of  Integration  Points 

100 

100 

100 

As  was  discussed  previously  the  FMM  and  MMM  are  very  similar,  the  main  differ¬ 
ence  being  that  the  weighting  functions  are  different.  Recall  the  FMM  uses  the  weighting 
functions  —  c’"",  with  v  being  the  angular  variable  in  the  new  coordinate  system  that 
was  generated  from  the  conformal  mapping.  An  interesting  question  is  whether  this 
choice  of  functions  is  better  than  the  set  of  functions  xpn  =  where  6  is  the  an¬ 
gular  variable  in  circular  coordinates.  If  this  were  true  the  conformal  mapping  would 
be  unnecesary,  since  the  entire  formulation  could  be  performed  in  circular  coordinates. 
However,  in  practice  the  FMM  requires  far  fewer  terms  than  by  using  V’n  =  c*”®.  For  ex¬ 
ample,  Figure  (4.1)  shows  the  backseat  ter  from  a  rigid  10:1  aspect  ratio  elliptic  cylinder 
for  rigid  boundary  conditions.  The  solid  curve  shows  the  result  using  the  FMM  and  16 
modes,  and  converges  to  the  exact  solution,  as  will  be  shown  in  the  following  section. 
The  dashed  curve  uses  the  cylindrical  weighting  functions  for  the  same  number  of  modes. 
The  dash-dotted  curve  used  the  cylindrical  weighting  functions  using  22  modes.  As  can 
be  seen  from  this,  the  cylindrical  weighting  functions  may  be  used,  but  will  require  more 
modes  than  the  FMM. 

4.2  Ellipse 

The  elliptic  cylinder  is  a  very  important  check  on  the  FMM  because  the  exact 
solution  is  know  in  terms  of  Mathieu  functions  (Appendix  B).  For  this  geometry,  the 
results  from  the  FMM  may  be  compared  with  the  exact  solution  for  both  impenetrable 
and  penetrable  fluid  cylinders  for  a  wide  range  of  aspect  ratios.  Results  will  be  presented 
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for  a  rigid  10:1  aispect  ratio  ellipse,  and  for  a  2:1  aspect  ratio  cylinder  for  both  rigid 
and  penetrable  boundary  conditions.  For  all  of  these  cases,  the  exact  solution,  auid  the 
solution  from  the  FMM  are  essentially  indistinguishable.  Figure  (4.2)  shows  the  geometry 
of  the  elliptic  cylinder  for  a  2:1  aspect  ratio. 

Figure  (4.3)  illustrates  the  backscatter  from  an  approximately  10:1  aspect  ratio 
rigid  elliptic  cylinder  vs.  kL  for  a  wave  incident  along  the  major  axis.  The  parameters 
used  are  shown  in  the  following  table: 


Mapping  Coefficients 

9i  = 

Aspect  Ratio 

10.033:1 

Number  of  Modes 

M=18 

Number  of  Integration  Points 

120 

Error 

CM  <  .0001 

Here  both  the  results  from  the  FMM  and  the  exact  result  from  separation  of  vari¬ 
ables  are  plotted  versus  kL,  As  can  be  seen,  the  results  are  indistinguishable.  Figure 
(4.4)  illustrates  both  the  exact  result  and  the  results  from  the  FMM  for  a  10:1  aispect 
ratio  ellipse  under  rigid  boundary  conditions  for  a  plane  wave  incident  along  the  minor 
axis.  Again,  the  results  are  indistinguishable. 

Figure(4.5)  illustrates  the  backscatter  from  an  approximately  2:1  aspect  ratio  rigid 
elliptic  cylinder  vs.  wavenumber.  The  parameters  used  were: 

Mapping  Coefficients  ?-i  =  l  =  e"’ 

Aspect  Ratio  2.164:1 

Number  of  Modes  M=12 

Number  of  Integration  Points  100 
Error  <  .0001 

Both  the  results  from  the  FMM  and  the  exact  result  from  separation  of  variables 
are  shown.  Again,  the  results  are  indistinguishable. 
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Figure  (4.6)  illustrates  the  backscatter  from  an  approximately  2:1  aspect  ratio 
penetrable  fluid  elliptic  cylinder  vs.  wavenumber.  The  material  parameters  used  were 
p! P\  =  2,  and  c/cj  =  1.05.  Here  both  the  results  from  the  FMM  and  the  exact  result 
from  separation  of  vau'iables  are  plotted.  Again,  the  results  are  indistinguishable,  hence 
illustrating  the  success  of  the  method. 

4.3  Square 

For  the  square  cylinder,  a  conformal  mapping  will  not  be  valid  at  the  sharp  corners. 
However,  after  truncating  the  mapping  the  corners  become  slightly  rounded,  and  the 
mapping  will  be  valid.  Therefore,  only  the  first  few  terms  in  the  mapping  will  be  used. 
The  parameters  used  for  the  square  cylinder  are: 

Mapping  Coefficients  5-i  =  l  53  =  —1/6  gj  =  1/56 

Number  of  Modes  M=30 

Number  of  Integration  Points  120 
Error  <  .0001 

The  cross  section  of  the  now  rounded  square  cylinder  appears  as  in  Fig.  (1.7). 
Results  for  the  physical  optics  approximation  as  well  as  the  FMM  will  be  generated  for 
the  rounded  square. 

Figure  (4.8)  illustrates  the  backscattered  form  function  for  a  plane  wave  incident 
at  45  degrees  from  a  rigid  square  cylinder.  The  backscattered  form  function  is  plotted 
versus  A'L,  where  L  is  the  length  of  a  side  of  the  square.  The  solid  curve  is  from  the  FMM, 
and  the  dashed  curve  is  from  the  physical  optics  method.  Note  that  at  low  frequency 
the  two  curves  do  not  coincide.  This  discrepancy  is  not  alarming  because  physical  ojitics 
is  valid  only  at  high  frequency.  As  the  frequency  increases,  the  agreement  between  the 
FMM  and  physical  optics  improves. 
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Similar  results  axe  obtained  for  a  plane  wave  incident  at  45  degrees  from  a  soft 
square  cylinder,  as  shown  in  Fig.  (4.9). 

Figure  (4.10)  illustrates  the  backscatter  of  a  plane  wave  incident  at  0  degrees  from 
a  rigid  square  cylinder.  Again,  the  two  curves  show  good  agreement  except  at  low 
frequency,  where  physical  optics  is  not  valid. 

For  the  penetrable  square  cylinder,  there  is  no  exact  solution.  However,  if  the 
scatterer  is  weak,  that  is  the  material  properties  of  the  scatterer  are  a  perturbation 
of  the  surrounding  fluid,  the  the  Distorted  Wave  Born  Approximation  (DWBA)  may 
be  used.  For  a  square  cylinder  with  ^  =  1.05,  ^  =  1.05  and  a  plane  wave  incident 
at  0  degrees,  the  scattering  integral  reduces  to  a  form  that  may  be  interpreted  as  the 
interference  between  a  reflection  off  the  front  surface  of  the  square  and  a  wave  that  is 
transmitted  through  the  front  face  of  the  cylinder,  reflects  off  the  back  face,  and  then  is 
transmitted  through  the  front  face.  By  summing  the  contributions  from  these  two  waves, 
the  DWBA  is  calculated  and  plotted  along  with  results  from  the  FMM  in  Fig.  (4.11). 
The  solid  curve  is  the  FMM  and  the  dashed  curve  is  the  DWBA.  Both  show  the  dramatic 
oscillatory  effect  due  to  the  interference  between  the  two  waves.  The  agreement  is  again 
excellent. 


4.4  Three  Leaf  Clover 

This  section  shows  a  more  complex  shape  than  considered  previously.  This  work  was 
conducted  because  it  is  fundamentally  different  than  the  ellipse  and  square  in  that  there 
are  concavities  in  the  surface.  The  geometry  is  shown  in  Fig.  (4.12).  The  parameters 
used  for  this  shape  are: 
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Figure  4.10:  Backscattered  form  function  of  a  Rigid  Square  :  =  0  (direction  of  in¬ 

cidence  from  the  side  of  the  square).  The  dashed  curve  is  is  from  the  physical  optics 
approximation,  and  the  solid  curve  is  from  the  FMM  under  soft  boundary  conditions. 
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Figure  4.11:  Backscattered  form  function  of  a  Weakly  Scattering  Square  :  =  0  (direc¬ 

tion  of  incidence  from  the  side  of  the  square).  The  dashed  curve  is  from  the  DWBA,  and 
the  solid  is  from  the  FMM.  The  material  parameters  are  ^  =  1.05  and  ^  =  1.05. 
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Mapping  Coefficients 

it-i— 1  §2  —  —-2  Qh  —  —.05 

Number  of  Modes 

M=30 

Number  of  Integration  Points 

120 

Error 

tM  <  .0001 

Figure  (4.13)  illustrates  the  backscatter  from  a  rigid  clover  for  a  plane  wave  incident 
at  ^0  =  ’'■/h  compared  with  physical  optics.  The  backscattered  form  function  is  plotted 
versus  kL,  where  L  is  the  average  radius  of  the  surface.  The  dashed  curve  is  physical 
optics,  and  the  solid  curve  is  from  the  FMM  under  rigid  boundary  conditions.  Again, 
physical  optics  is  not  valid  at  low  frequencies,  so  the  comparison  there  is  not  expected  to 
be  close,  but  at  high  frequencies  the  comparison  improves.  For  both  this  case,  and  that 
of  the  square  cylinder,  even  at  high  frequencies  the  agreement  is  close  but  not  perfect. 
This  is  to  be  expected  because  physical  optics  assumes  that  each  point  scatters  as  though 
it  were  part  of  an  infinite  plane,  which  is  not  true.  It  also  ignores  any  contribution  from 
circumferential  waves,  which  tend  to  alter  the  structure  of  the  scatter  pattern  .  Figure 
(4.14)  shows  physical  optics  (solid  curve)  compared  with  the  exact  solution  for  a  2:1 
aspect  ratio  elliptic  cylinder.  The  dashed  curve  is  for  soft  boundary  condtions,  and  the 
dotted  curve  is  for  rigid  boundary  conditions.  As  can  be  seen,  the  physical  optics  method 
will  not  necesarily  accurately  predict  the  level  of  the  nulls. 
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kL 


Figure  4.14:  Backscattered  form  function  from  a  2:1  Ellipse.  Oo  =  0.  The  solid  curve  is 
is  from  the  physical  optics  approximation,  the  dashed  curve  is  from  the  exact  solution 
under  soft  boundary  conditions,  and  the  dotted  curve  is  from  the  exact  solution  under 
rigid  boundary  conditions.  This  illustrates  the  fact  that  the  phj'sical  optics  solution  will 
not  necesarily  accurately  predict  the  level  of  the  nulls. 
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Chapter  5 


Discussion  and  Conclusions 


5.1  Conclusions 

The  development  of  a  new  method  for  calculating  the  farfield  scattered  pressure  of  a 
plane  wave  due  to  a  cylinder  of  noncircular  cross  section  has  been  presented  in  this  thesis. 
The  new  method  has  been  shown  to  be  very  accurate  and  numerically  efficient  for  a  wide 
range  of  cross  sections  and  boundary  conditions.  The  formulation  is  simple  in  that  it  does 
not  need  to  be  tailored  to  the  shape  of  the  scatterer,  as  does  the  T-Matrix.  Numerical 
results  presented  for  elongated  cross  sections,  as  well  as  cross  sections  with  sharp  edges 
and  concavities,  have  shown  excellent  agreement  with  exact  and  approximate  solutions. 
It  has  also  been  shown  to  be  numerically  more  efficient  than  competing  approaches  when 
the  scatterer  is  not  close  to  circular. 

This  method  most  nearly  resembles  the  MMM,  in  that  the  scattered  field  is  de¬ 
scribed  as  a  finite  summation  of  nonorthogonal  modal  functions.  The  difference  is  in 
how  the  scattered  field  coefficients  are  chosen.  In  the  MMM,  the  mean  squared  error  on 
the  boundary  is  minimized.  The  drawback  of  this  approach  is  that  in  order  to  satisfy  the 
boundary  conditions,  far  more  coefficients  are  required  than  are  necesary  to  describe  the 
farfield  scattered  pressure. 
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It  is  here  that  the  method  developed  in  this  thesis,  the  FMM,  becomes  useful. 
Instead  of  minimizing  the  mean  squared  error,  the  scattered  field  coefficients  are  chosen 
so  that  there  is  no  error  in  the  Fourier  components  e*""  for  |n|  <  N.  In  practice  fewer 
terms  are  required  to  have  a  convergent  farfield  scattered  pressure  than  to  satisfy  the 
boundary  conditions  rigorously.  Therefore,  if  nearfield  quantities  are  not  required,  the 
FMM  will  in  general  be  more  eflBcient  than  the  MMM. 

Another  benefit  of  this  research  is  the  fact  that  the  conformal  mapping  portion  may 
be  used  in  conjunction  with  the  T-Matrix  to  improve  the  range  of  surfaces  the  method  can 
handle.  For  example,  by  expressing  the  surface  fields  in  terms  of  the  new  angular  variable 
y,  accurate  results  were  achieved  for  an  elongated  (elliptical)  cross  section  using  the  T- 
Matrix,  where  the  Green’s  function  was  expressed  in  terms  of  cylindrical  eigenfunctions. 
This  improvement  was  due  to  the  fact  that  sampling  points  for  the  numerical  integrations 
were  more  rapidly  sampled  at  the  ends  of  the  ellipse,  therefore  eliminating  the  error  that 
causes  the  method  to  become  unstable. 

Although  the  FMM  may  be  more  efficient  and  general  than  both  the  MMM  and 
T-Matrix,  the  major  disadvantage  of  the  FMM  is  that  it  has  not  been  rigorously  shown  to 
converge  to  the  true  solution.  Results  have  been  presented  that  show  excellent  agreement 
with  reference  solutions.  By  proof  of  example.  The  author  is  confident  that  the  solution 
is, in  fact,  correct.  In  order  to  use  the  FMM,  it  is  necesary  to  change  variables  to  a  new 
coordinate  system  in  which  the  radial  coordinate  u  =  0  is  the  scatterer  surface.  This  is 
achieved  through  a  conformal  mapping,  namely  conformally  mapping  the  exterior  of  a 
closed  surface  to  a  semi-infinite  strip,  and  the  interior  of  a  closed  surface  to  a  semi-infinite 
strip.  To  the  knowledge  of  the  author,  this  is  the  first  time  that  an  efficient  algorithm 
has  been  presented  for  doing  this. 

In  conclusion,  a  new  technique  for  solving  the  problem  of  scattering  from  a  cylinder 
of  noncircular  cross  section  has  been  presented.  This  technique  has  proven  to  be  accurate 
for  a  wide  variety  of  cross  sections  and  boundary  conditions,  as  well  as  numerically 
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eflBcient.  Severai  examples  of  baokscatter  calculations  were  presented  and  compared  to 
either  exact  or  approximate  solutions.  As  discussed  in  the  following  section,  however, 
some  interesting  research  remains  to  be  done. 


5.2  Recommendations  for  Future  Research  ri 

There  are  several  types  of  research  that  could  follow  directly  from  the  work  pre¬ 
sented  in  this  thesis.  Since  an  efficient  algorithm  has  been  developed  to  describe  the 
scattering  from  a  noncircular  cylinder,  it  would  be  straightforward  to  generate  time- 
domain  scattering  data  from  which  the  physics  underlying  the  scattering  process  could 
be  studied. 

Although  the  focus  of  this  research  has  been  the  scattering  of  sound  by  a  (pene¬ 
trable)  fluid  cylinder,  the  extension  to  an  elastic  cylinder  would  be  straightforward.  For 
example,  the  shear  wave  that  may  be  supported  in  an  elastic  scatterer  may  be  accounted 
for  by  introducing  another  set  of  unknown  coefficients  for  that  new  field,  and  then  that 
set  of  coefficients  would  be  eliminated  by  introducing  the  extra  boundary  condition  that 
arises  for  an  elastic  boundary. 

A  third  area  that  would  be  interesting  to  pursue  would  be  that  of  a  three-dimensional 
scatterer.  This  would  be  more  difficult  because  conformal  mapping  involves  only  two  di¬ 
mensions  and  might  not  be  extendable  to  three  dimensions. 

A  fourth  area,  alluded  to  earlier  in  the  thesis,  involves  using  the  fact  that  the 
Helmholtz  equation  in  the  new  coordinates  is  the  same  as  the  Helmholtz  equation  in 
Cartesian  coordinates  with  a  spatially  dependent  wavenumber.  Propagation  of  an  acous¬ 
tic  wave  in  a  range  and  depth  dependent  environment  is  an  important  and  difficult 
problem,  and  the  methods  developed  in  this  thesis  may  give  insight  into  a  new  technique 
for  solving  this  problem. 
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Chapter  6 


Appendices 


6.1  Appendix  A:  Numerical  Solution  of  Nonlinear  System  of 
Equations 

This  appendix  describes  how  to  numerically  solve  the  set  of  nonlinear  equations 
that  arise  from  the  conformal  mapping  procedure.  The  procedure  is  a  simple  extension 
of  the  Newton-Raphson  procedure  [35].  The  Newton- Raphson  method  is  an  iterative 
algorithm  that  linearizes  the  nonlinear  system  of  equations.  For  example,  consider  the 
one-dimensional  case 

d(6)  =  0  (6.1) 

where  d{6)  is  a  nonlinear  function  of  6,  and  6  is  the  unknown.  Assume  an  estimate  of 
the  value  of  the  root  of  the  equation  is  6®,  such  that  d{6^  -f  /i®)  =  0,  where  /i®  is  small. 
This  may  be  written  as  a  Taylor  series 

+  ft®)  =  d(ftO)  +  ft®d'(6®)  +  •  •  •  (6.2) 

Retaining  only  the  first  term,  and  setting  d(6°  -}-  ft®)  =  0,  the  correction  term  may  be 
solved  for;  ft®  =  d{S^)/d'{6°).  The  new  estimate  of  the  root  thus  becomes  =  6®  -(-  /»®. 
The  entire  procedure  is  now  repeated,  with  the  Mth  correction  term  being  given  b} 
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*  =  d{6^  ^)ld'{8^  *).  In  practice  M  may  be  chosen  such  that  d{6*^)  is  less  than 


some  tolerance. 


This  is  easily  generalizable  to  the  N-dimensional  case.  The  problem  in  N-dimensions 


>  ^n)  —  0 


-I^n)  =  0  (6.3) 

Following  the  same  procedure  as  in  the  one-dimensional  case,  the  Mth  correction  vector 
is  obtained  from  solving  the  following: 

=  0  (6.4) 

where  is  the  A/th  correction  vector,  and  is  the  Jacobian  matrix  whose  ele¬ 


ments  are  given  by 


jM-l  _  I 


Solving  for  ' ,  a  new  estimate  for  the  solution  vector  becomes: 


-I- 


Again,  this  process  is  iterated  until  the  maximum  of  d(6^^)  is  less  than  some  tolerance. 

Note  here  that  most  of  the  computation  time  in  the  Newton- Raphson  method 
is  in  calculating  the  Jacobian  matrix.  The  modified  Newton-Raphson  avoids  this  by 
calculating  the  Jacobian  only  for  the  original  estimate  of  the  solution  vector  6^  and 
solving  the  following  equation 

d(^^’^)  + JV'-' =  0  (6.7) 


Using  this  method  will  generally  require  a  larger  M.  but  it  is  much  faster  because  the 
Jacobian  only  has  to  be  calculated  once. 


The  convergence  in  either  method  will  be  very  fast  provided  the  initial  guess  ^ 
is  close  to  the  true  solution.  For  shapes  that  are  nearly  circular,  a  reasonable  initial 
guess  is  that  ^  =  Q.  For  shapes  that  differ  considerably  from  circular,  an  extension  of 
Newton- Raphson  must  be  tried.  Recall  that  on  the  surface 

»  +  +  E  (6.8) 

n=l 

If  the  Rn  are  not  small,  the  new  procedure  is  to  divide  all  Rn  by  some  suitably  large 
integer  L,  leading  to  the  sequence  of  problems 

X  +  =  ae'*  +  I  f;  +  fine"'-'*]  (6.9) 

The  problem  for  /  =  1  is  solved  using  the  intial  guess  of  ^  =  0.  Once  the  ^  have  been 
solved  for.  these  values  are  used  as  the  initial  guess  for  I  —  0,.  The  procedure  is  then 
continued  until  I  =  X,  which  is  the  original  problem. 
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6.2  Appendix  B:  Equations  used  in  Numerical  comparisons 


This  appendix  will  describe  the  equations  used  in  generating  the  numerical  results 
various  existing  solutions.  The  results  are  used  in  the  main  text  for  comparison  with  the 
FMM.  Included  are  descriptions  of  the  1)  exact  solutions  which  involves  separation  of 
variables  for  the  two  surfaces  which  are  applicable  to  this  thesis,  circular  and  elliptical, 
for  both  penetrable  and  impenetrable  scatterers,  2)  the  physical  optics  approximation  for 
high  frequency  impenetrable  scatterers,  and  3)  the  Distorted  Wave  Born  Approximation 
for  weakly  scattering  (highly  penetrable)  bodies  for  all  frequencies. 

6.2.1  Exact  Solutions  for  Separable  Geometries 


As  discussed  earlier,  the  conformal  mapping  changes  the  Helmholtz  equation  in 
Cartesian  coordinates  to  a  new  Helmholtz  equation 

F’uu  +  +  Ar^^(t£,u)jP  =  0  (6.1) 

where  /■(u,u)  depends  on  the  specific  mapping.  If  J^{u,v)  is  a  sum  of  a  function  of  u 
and  a  function  of  u,  then  the  equation  is  said  to  be  separable.  That  means  that  it  may 
be  solved  by  separation  of  variables.  The  equation  becomes: 

F’uu  +  Pvv  +  [/(«)  +  T*  =  0  (6.2) 


By  assuming  that  the  solution  is  a  product  of  a  function  of  u  and  a  function  of  r. 
P  =  U{u)V{v)^  and  insertng  this  into  Eq.(6.2),  leaves 

^  +  A:V(«)  =  “  (6.3) 


Since  the  left  side  of  this  equation  is  a  function  of  xi  only,  and  the  right  side  is  a  function 
of  V  only,  each  side  must  be  constant.  This  is  known  as  a  separation  constant  and  is 
denoted  by  7.  The  equation  now  becomes 


U 


+  k\f{u)  =  - 


=  7 


(6.4) 
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This  may  now  be  rewritten  as  two  equations 


V^^  +  k^givW  +  'rV  =  0 

Uuu  +  =  0  (6.5) 

These  are  ordinary  differential  equations  and  thus  may  be  solved  rigorously.  In  general 
the  first  of  these  equations  has  a  distinct  solution  for  infinitely  many  values  of  7.  Each 
of  these  values  7^  is  called  a  characteristic  value  or  eigenvalue,  and  the  corresponding 
solution  I4i(t;)  is  called  an  eigenfunction.  Once  the  eigenvalue  has  been  determined,  the 
second  equation  may  be  solved  to  yield  f4„(u).  The  total  solution  is  then  a  summation 
of  each  of  these  solutions,  i.e. 

CO 

£  h^V^{v)U„,{u)  (6.6) 

m=— 00 

Circular 

The  simplest  separable  geometry  is  a  circular  boundary.  For  this  case,  /(u)  =  e^“,  and 
g{v)  =  0.  The  resulting  ordinary  differential  equations  become 

Kv  +  7V  =  0 

-  7C/  =  0  (6.7) 

Solving  these  differential  equations  yields  the  following  solutions 

00  _ 

P"‘=  E  +  (6.8) 

m=— oo 

Where  the  Om  are  the  known  incident  coefficients,  and  are  the  unknown  scattered  field 
coefficients.  In  the  interior  of  the  scatterer,  the  corresponding  solution  is 

OO 

p-‘=  ^  d,„e‘’"V„(A7e“)  (6.9) 

Tn=— OO 

Where  the  dm  are  unkown  transmitted  field  coefficients.  The  boundary  conditions  are 

^ezt  _  ^int 

=  n-VP'”'  (6.10) 
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Inserting  Eqs.(6.8  and  6.9)  into  Eqs.(6.10),  and  utilizing  the  orthogonality  properties  of 
the  angular  eigenfunctions  leaves 

(6.11) 

Pi  Cl 

This  is  precisely  the  solution  obtained  from  the  formulation  in  Chapter  2. 

Elliptic 

The  next  simplest  separable  geometry  is  the  elliptic  boundary.  For  this  case,  /(u)  = 
2cosh(2u),  and  g(v)  =  — 2cos(2«).  The  resulting  ordinary  differential  equations  become 

-  2k^  cos(2v)V  +  =  0 

Uuu  +  cosh(2u)f/  —  yU  =  0  (6.12) 

The  first  of  these  equations  is  Mathieu’s  equation,  and  the  second  is  Mathieu’s  modified 
equation.  Note  that  Mathieu’s  equation  depends  on  wavenumber,  which  is  fundamentally 
different  than  the  corresponding  circular  case,  where  the  angular  equation  was  indepen¬ 
dent  of  wavenumber.  This  makes  no  difference  for  the  impenetrable  case,  but  for  the 
penetrable  scatterer  this  means  that  there  are  no  orthogonality  relations  between  the 
exterior  and  interior  angular  eigenfunctions. 

Solving  these  differential  equations  yields  the  following  solutions  in  McLachlan’s 
notation  [10]: 


OO 


^  [^2m^^2m{k,u)ce2„,{k,v)  +  al^^^Se2,n+2{k,u)se2m+2{k,v)  -f 

m=0 

^)c62Tn+l(^  *!  ^0  ^)] 


OO 


+  12  [blm^4mik,u)ce2^{k.  v)  +  6L+2-^^4m+2(^^»)-^C2m+2(^,v)  + 

m=:0 

b2m+i^^2m+lik,u)ce2m+l{k..  r)  -f-  6L+1^4m+l  l^)-^C2m  +  l  v)] 


(6.13) 
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Where  the  a  are  the  known  incident  coefficients,  and  the  6  are  the  unknown  scattered 
field  coefficients.  The  superscripts  and  are  for  even  or  odd  solutions,  and  the 
subscripts  2m,  2m+i  tcH  whether  the  angular  function  is  periodic  in  x  or  in  2x.  Therefore 
there  are  4  distinct  types  of  Mathieu  functions. 


1.  Even  in  2ir  : 

2.  Odd  in  2x  :  se2m+i(^,f) 

3.  Even  in  x  :  ce2m(^,v) 

4.  Odd  in  x  :  se2m(^. u) 


The  ce{k,v),se{k,v)  are  angular  Mathieu  functions  corresponding  to  cosine  and 
sine  functions,  respectively.  The  Me^^\k,u),  Ne^^\k,u)  are  radial  Mathieu  functions 
corresponding  to  Hankel  functions  of  the  first  kind,  and  Ce{k,u),Se{k,u)  are  radial 
Mathieu  functions  corresponding  to  Bessel  functions.  The  interior  solution  is 

OO 

pxnt  _  ld2mCe2m(kl,u)ce2mif^l,v)  +  <^7,r,j.oSe2m+2(kl,u)se2mA-2(kj,  v)  + 
m=0 

^m+1  ^®2m+l (^1 1  ^)ce2m+l (^1 1  v)  -f  d2,„^j5'e2m+l (^1 1  ^)5e2m+l (^1 » 

Where  the  d  are  the  unknown  transmitted  field  coefficients. 


Since  Mathieu  functions  are  not  available  in  MATLAB,  they  had  to  be  written. 
They  were  coded  directly  from  McLachlan,  with  the  radial  Mathieu  functions  being 
expanded  in  terms  of  products  of  Bessel  functions. 


For  the  rigid  cylinder,  the  boundary  condition  is 


n  •  =  0 


Using  the  orthogonality  relations  leaves 


_  Ce'2r„{k,uo) 

2m  -  ,(i)  “2m 


(6.15) 
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Since  the  interior  wavenumber  is  now  different  than  the  exterior  wavenumber,  ther  is  no 
longer  orthogonality  between  the  interior  and  exterior  angul2ir  wavefunctions.  By  multi- 
plyng  by  the  exterior  angular  eigenfunctions  and  integrating,  the  result  is  a  summation 
on  the  right  side  of  the  equation: 

OO 

'  (^1  ^o)  “t"  ^o)  ~  >  ^o) 

n=0 

65„Me'‘»(euo)  +  4^Ce',„{k,  u„)  =  £  a,.„di„Ce',Jk,,u„)  (6.18) 

n=0 

r2x 

Oln,m  =  /  ce2n{k,uo)ce2m{ki,uo)  (6.19) 

Jo 

This  is  equivalent  to  expanding  each  ce2m{ki,uo)  in  terms  of  ce2m{k,uo)-,  which  is  what 
Yeh[ll]  did.  This  is  only  for  even  solutions  periodic  in  tt,  and  there  are  corresponding 
results  for  the  other  types  of  solutions.  These  equations  may  be  solved  for  62,,,,  but  now 
each  is  coupled  to  all  of  the  a|„,,  i.e. 

00 

=  (6-20) 

n=0 

Similar  expressions  may  be  derived  for  the  remaining  scattered  field  coefficients. 


SO 


6.2.2  High  Frequency  (Physical  Optics)  Approximation 

The  backscattered  form  function  for  the  physical  optics  approximation  is  given  by 
the  following  integral  [19] 

rM  =  l/k  ■  (6.21) 

where  is  a  unit  vector  in  the  direction  of  the  incident  wave,  n  is  the  unit  normal  vector, 
k  is  the  incident  wavevector,  f  is  the  position  vector  of  each  point  on  the  surface,  dA  is 
a  differential  surface  element,  and  A  is  the  ensonified  area. 

Because  the  surface  normal  n  is  very  easy  to  calculate  in  the  new  coordinate  system, 
the  integral  will  be  evaluated  in  this  coordinate  system.  Using  the  properties  of  the  new 
coordinate  system,  the  quantities  in  the  integral  become; 

.  n  =  ■i[i„cos(^o)  +  l/u  sin(^o)] 

*l%i, 

k  •  r=  k{x  cos(0o)  +  y  sin(^o) 

dA=  Kdv  (6.22) 


Inserting  these  into  the  above  integral  leaves 


All  of  the  quantities  appearing  in  the  above  integral  are  easily  calculated  once  the  map¬ 
ping  coefficients  are  known.  For  example,  x  and  y  are  found  by  taking  the  real  and 
imaginary  parts  of  the  mapping  function: 


xAiy  =  g-^c"’  (6-24) 

TinO 
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a:„  and  may  be  found  by  taking  the  real  and  imaginary  parts  of  the  derivative  of  the 
mapping  function  with  respect  to  u: 

OO 

Xu  +  iyu  =  g-ie""  -  (6.25) 

n=0 

The  mapping  coefficients  used  for  the  shapes  in  this  thesis  have  already  been  de¬ 
termined  elsewhere  in  this  thesis,  and  are  summarized  as  follows: 

•  10:1  ellipse  :  g-i  =  1,  =  e~-^ 

•  2:1  ellipse  :  g^i  =  1, 

•  square  :  g_i  =  1,  ^3  =  -1,  ^7  =  ^ 

•  clover  :  g_i  =  1,  ^2  =  --2,  gs  =  -  05 

6.2.3  Distorted  Wave  Born  Approximation 

For  weakly  scattering  objects,  that  is  scatterers  whose  material  properties  are  a 
perturbation  of  the  surrounding  fluid,  the  Born  approximation  assumes  the  total  field 
on  the  scatterer  surface  to  be  the  incident  field,  and  then  uses  Green’s  theorem  to  gen¬ 
erate  the  scattered  pressure  by  integrating  over  the  surface.  An  extension  to  the  Born 
approximation  found  using  the  T-Matrix  is  given  by  Waterman  [27].  His  result  is  for 
three-dimensional  scatterers.  but  for  two  dimensions  it  is  given  by: 

F-“‘(7r)  [(p,  -  p)  +  (c,  -  c)]  ^  (6.26) 

where  dS  is  differential  surface  element.  For  the  square  cylinder  and  plane  wave  incident 
along  the  side  of  the  square,  performing  this  integration  yields 

[1  -  (6.27) 
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The  physical  interpretation  of  this  result  is  as  follows:  the  first  term  in  the  brackets 
represents  the  reflection  of  the  incident  wave  off  the  front  face.  The  second  term  represents 
a  wave  that  is  transmitted  through  the  front  face  of  the  scatterer,  is  reflected  oflT  the  back 
face,  and  then  is  transmitted  back  through  the  front  face.  This  can  be  seen  from  the 
fact  that  the  phase  of  the  transmitted  wave  is  k2L,  and  the  length  of  the  path  described 
above  is  2L.  However,  the  wavenumber  in  the  scatterer  is  not  k,  but  ki .  The  Distorted 
Wave  Born  Approximation  is  found  by  accounting  for  this  fact.  That  is,  the  wave  that  is 
transmitted  through  the  cylinder  will  have  a  wavenumber  fci,  where  k\  =  kf^.  Replacing 
k  by  k^  in  the  transmitted  term  yields  for  the  farfleld  pressure: 

/>'“■(!)  !((.,  -  P)  4  (c,  -  c)|  L,'*'-  (l  -  (6.28) 

In  order  to  apply  this  method  to  a  general  scatterer,  it  is  necesary  to  be  able  to 
evaluate  the  original  integral  in  closed  form,  and  then  to  separate  out  the  contributions 
of  the  waves  that  are  propagating  on  the  interior  of  the  scatterer.  In  general  this  is  not 
possible,  so  this  method  can  only  be  applied  in  certain  simple  cases. 
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